Renormalization vs Strong Form Factors for One Boson Exchange Potentials. 
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We analyze the One Boson Exchange Potential from the point of view of Renormalization theory. 
We show that the nucleon-meson Lagrangean while predicting the NN force does not predict the 
NN scattering matrix nor the deuteron properties unambiguously due to the appearance of short 
distance singularities. While the problem has traditionally been circumvented by introducing vertex 
functions via phenomenological strong form factors, we propose to impose physical renormalization 
conditions on the scattering amplitude at low energies. Working in the large N c approximation with 
n,a,p and ui mesons we show that, once these conditions are applied, results for low energy phases 
of proton-neutron scattering as well as deuteron properties become largely insensitive to the form 
factors and to the vector mesons yielding reasonable agreement with the data and for realistic values 
of the coupling constants. 
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I. INTRODUCTION 



The One Boson Exchange (OBE) Potential has been 
a cornerstone for Nuclear Physics during many years. It 
represents the natural generalization of the One Pion Ex- 
change (OPE) potential proposed by Yukawa [l[, and also 
the scalar meson potential introduced by Johnson and 
Teller [2| . With the advent of vector mesons these degrees 
of freedom were included as well BELIE H- Actually, 
Regge theory yields such a potential within a suitable 
approximation Q. The disturbing short distance diver- 
gences were first treated by using a hard core boundary 
condition [1, 0, Q and it was soon realized that diver- 
gences in the potential could be treated by introducing 
phenomenological form factors mimicking the finite nu- 
cleon size [1. The field theoretical OBE model of the NN 
interaction^, [l(| includes all mesons with masses below 
the nucleon mass, i.e., n, rj, p(770) and a; (782). We refer 
to Ref. [ill, EH for accounts of the many historical re- 
volves of the problem. An important lesson from these 
developments has been that the non-perturbative nature 
of the NN force is better handled in terms of quantum me- 
chanical potentials at low energies where relativistic and 
nonlocal effects contribute at the few percent level. Al- 
though such a framework has remained a useful, appeal- 
ing and accurate phenomenological model after a suit- 
able introduction of phenomenological strong form fac- 
tors [13, EH it is far from being a complete description of 
the intricacies of the nuclear force. The highly successful 
Partial Wave Analysis (PWA) of the Nijmegen group 14 1 
while providing a spectacular fit with % 2 /DOF < 1 com- 
prising a large body of pn and pp scattering data checks 
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mainly OPE and some contributions from other mesons, 
since the interaction below 1.4fm is parameterized by an 
energy dependent square well potential. 

A traditional test to NN forces in general and OBE 
potentials in particular has been NN scattering in the 
elastic region. In such a situation relative NN de Broglie 
wavelengths larger than half a fermi are probed; a factor 
of two bigger scale than the Compton wavelengths of the 
vector and heavier mesons. However, while from this sim- 
ple minded argument we might expect those mesons to 
play a marginal role, OBE potentials have traditionally 
been sensitive to short distances requiring a unnatural 
fine tuning of the vector meson coupling. As a conse- 
quence there has been some inconsistency between the 
couplings required from meson physics, SU(3) or chiral 
symmetry on the one hand and those from NN scatter- 
ing fits on the other hand (see also [IE, El, E3)- Part 
of the disagreement could only be overcome after even 
shorter scales were explicitly considered [H, EH . A more 
serious shortcoming stems from the use of strong form 
factors which have mainly been phenomenologically mo- 
tivated and loosely related to the field theoretical meson- 
baryon Lagrangean from which the meson exchange pic- 
ture is derived. It is therefore not exaggerated to say 
that strictly speaking the OBE potentials have not been 
solved yet. Of course, this may appear as a mathemati- 
cally interesting problem with no relevance to the physics 
of NN interactions. However, as we will see, the Meson- 
Nucleon Lagrangean itself while providing the NN OBE 
potential from the Born approximation, does not predict 
the NN S-matrix and the deuteron unambiguously be- 
yond perturbation theory from the OBE potential. The 
unspecified information in the Lagrangean can be ad- 
vantageously tailored to fit the data in the low energy 
region. We will also show that once this is done, the ver- 
tex functions play a minor role, with a fairly satisfactory 
description of central waves and the deuteron. 
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It is notorious that the OBE potentials although expo- 
nentially suppressed with the corresponding meson mass, 
~ e~ mr , are by themselves large at short distances and 
mostly even diverge as 1/r 3 . For a singular potential, 
i.e. a potential fulfilling lim r ^o 2/i\V(r)\r 2 = oo [20. l2l|. 
the Hamiltonian is unbounded from below preventing 
the existence of a stable two nucleon bound state when 
lim r _>o 2fiV(r)r 2 < —1/4. Of course, the singularity is 
unphysical as it corresponds to the interaction of two 
point-like static classical particles and not to extended 
nucleons with a finite size of about half a fermi 1 . From 
a quantum mechanical viewpoint, however, the relative 
NN de Broglie wavelength provides the limiting resolu- 
tion scale physically operating in the problem. This sug- 
gests that finite nucleon size effects should also play a 
marginal role in NN scattering in which case it should be 
possible to formulate the problem without any explicit 
reference to form factors. In fact, as we will explicitly 
demonstrate Renormalization is the natural mathemati- 
cal tool to implement the physically desirable decoupling 
of short distance components of the interaction at the 
energies involved in NN elastic scattering. 

Within the NN system the problem of infinities has 
traditionally been cured [3] by the introduction of phe- 
nomenologically or theoretically motivated strong form 
factors in each meson- nucleon vertex, T m ^N(q 2 ) (m — 
ir, a, p, uj etc.) where the off-shellness of the nucleon legs 
is usually neglected. This procedure somewhat mimics 
the finite nucleon size but strong form factors are fitted 
and constrained in practice to NN scattering data and 
deuteron properties. This corresponds to the replace- 
ment of the potential V m (q) — > V m (q) [r TO jv,/v(<7)] 2 where 
typically a monopole form is taken for each separate me- 
son T mNN (q) = (Al lNN - m 2 )/(A mNN - q 2 ) and gen- 
erally A m NN ~ 1 — 2GeV 2 . Due to the long distance 
distortion introduced by the vertex function deuteron 
properties impose limitations on the lowest cut-off value 
A^nn > 1.3GeV still fitting the result [13]. 

Because of their fundamental character and the crucial 
role played in NN calculations there have been count- 
less attempts to evaluate strong form factors by sev- 
eral means, mainly T^nnU^)- These include meson 
theory [H, H, H, M M J27| jRegge models [H, chi- 
ral soliton models [13, [13 EH, H3], QCD sum rules [33] 
the Goldberger-Treimann discrepancy [341 ] or lattice 
QCD [35|, H3] and quark models [37| ■ Most calculations 
yield rather small values A^nn ~ 800MeV generating 
the soft form factor puzzle for the OBE potential for sev- 
eral years since the cut-off could not be lowered below 
Attnn — 1.3GeV without destroying the quality of the 
fits and the description of the deuteron [13] ■ The contra- 



1 Of course, the nucleon size depends on the particular electroweak 
probe. We give here a typical number. 

2 For a monopole the operating scale is lower, A m jviv /%/2, because 
the square of the form factor enters in the modification of the 
potential. 



diction was solved by including either pir exchange [38| . 
a strongly coupled excited 7r'(1300) state [13], two pion 
exchange [13] or three pion exchange [!l|] . Some of these 
ways out of the paradox assume the meson exchange pic- 
ture seriously to extremely short distances. However, as 
noted in Ref. [13] the contradiction is misleading since 
a large cut-off is needed just to avoid a sizable distor- 
tion of the OBE potential in the region r > 0.5fm which 
can also be achieved by choosing a suitable shape of the 
form factor. This point was explicitly illustrated by us- 
ing the Skyrme soliton model form factors [13]. In fact, 
this conclusion is coherent with the early hard core regu- 
larizations 0, 0, , recent lattice calculations [13] (where 
an extremely hard A^nn ~ 1.7GeV and a rather flat be- 
haviour are found) and, as we will show, with the renor- 
malization approach we advocate. 

The implementation of vertex meson-nucleon functions 
has also notorious side effects, in particular it affects 
gauge invariance, chiral symmetry and causality via dis- 
persion relations. As it is widely accepted, besides the 
description of NN scattering and the deuteron, one of 
the great successes and confirmations of Meson theory 
has been the prediction of Meson Exchange Currents 
(MEC's) for electroweak processes (see [13, EH for re- 
views and references therein) . In the case of gauge invari- 
ance, the inclusion of a form-factor introduced by hand, 
i.e., not computed consistently within meson theory, im- 
plies a kind of non-locality in the interaction. This can be 
made gauge invariant by introducing link operators be- 
tween two points, thereby generating a path dependence, 
and thus an ambiguity is introduced. In the limit of weak 
non-locality the ambiguity is just the standard operator 
ordering problem, for which no obvious resolution has 
been found yet. Form factors can also be in open conflict 
with dispersion relations, particularly if they imply that 
the interaction does not vanish as a power of the momen- 
tum everywhere in the complex plane. We will show that 
within the renormalization approach, all singularities fall 
on the real axes and spurious deeply bound states are 
shifted to the real negative. The extremely interesting is- 
sue of analyzing the consequences of renormalization for 
electroweak processes is postponed for future research. 

In the present paper we approach the NN problem for 
the OBE potential from a renormalization viewpoint. We 
analyze critically the role played by the customarily used 
phenomenological form factors. As a viable alternative 
we carry out the renormalization program to this OBE 
potential to manifestly implement short distance insensi- 
tivity as well as completeness of states by removing the 
cut-off. In practice, we use the coordinate space renor- 
malization by means of boundary conditions pit l45l |46| . 
The equivalence to momentum space renormalization us- 
ing counterterms for regular and singular potentials was 
discussed in Refs. [47l l48j] . In order to facilitate and 
simplify the analysis we will use large N c relations for 
meson-nucleon couplings [13, H3 HH which are well sat- 
isfied phenomenologically and pick the leading tensorial 
structures for the OBE potential. In this picture mesons 
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are stable with their mass scaling asm~ N®, nucleons 
are heavy with their mass scaling as ~ N c , the NN 
potential also scales as Vnn ~ N c . The OBE compo- 
nent is dominated by the n, a, p and w mesons IHiij. 
Further advantages of using this large N c approximation 
have been stressed in in regard to Wigner and Serber 
symmetries in Refs. [HJ, [H, [Hy, H3| , in particular the fact 
that relativistic, spin-orbit and meson widths corrections 
are suppressed by a relative factor suggesting a 

bold 10% accuracy. However, we hasten to emphasize 
that despite the use of this appealing and simplifying ap- 
proximation in the OBE potential we do not claim to 
undertake a complete large N c calculation since multiple 
meson exchanges and A intermediate states should also 
be implemented [53|. In spite of this, some of our re- 
sults fit naturally well within naive expectations of the 
large N c approach. The coordinate space renormaliza- 
tion scheme is not only convenient and much simpler but 
it is also particularly suited within the large N c frame- 
work, where non-localities in the potential are manifestly 
suppressed, and an internally consistent multimeson ex- 
change scheme is pos sible if energy independent poten- 
tials are used [58l.l59||. 

The paper is organized as follows. In Sect. [IT] we write 
down a chiral Lagrangean in order to visualize the calcu- 
lation of the OBE potential in the large N c limit and an- 
alyze its singularities. The standard approach to prevent 
the singularity has been to include form factors to rep- 
resent vertex functions, an issue which is analyzed crit- 
ically in Section IIIII where the alternative between fine 
tuning and the appearance of spurious bound states is 
highlighted. In Section HVl we discuss the physical condi- 
tions under which a description of NN scattering makes 
sense within a renormalization point of view. In addi- 
tion, we discuss some general features which apply to the 
solutions of the Schrodinger equation with the local and 
energy independent large iV c OBE potential on the basis 
of renormalization. In Sect. [V] we- analyze the So chan- 
nel from which the scalar meson parameters may be fixed. 
We also discuss the role played by spurious bound states 
which appear in this kind of calculations. The deuteron 
and the corresponding low energy parameters as well as 
the 3 Si — 3 D\ phase shifts are analyzed in Sect. [Vj The 
marginal influence of form factors in the renormalization 
process is shown in Section IVIII Finally, in Sect. IVIIII 
we summarize our main points and conclusions. In ap- 
pendix [A] we also review current values for the coupling 
constants from several sources entering the potential. 



II. OBE POTENTIALS AND THE NEED FOR 
RENORMALIZATION 

In this section we briefly sketch the well known process 
of deriving the OBE potential from the nucleon-meson 
Lagrangean. We appeal a chiral Lagrangean as done 
in Refs. [HI, [TH, [TtJ and keep only the leading N c con- 
tributions to the OBE potential due to the tremendous 



simplification which proves fair enough to illustrate our 
main point, namely the lack of uniqueness of the S- matrix 
from the OBE potential. In a more ellaborated version 
the present calculation should include many other effects 
such as relativistic corrections, spin-orbit coupling, me- 
son widths, multi-meson exchange and A intermediate 
states. 



A. Meson-Nucleon Chiral Lagrangian 

We use a relativistic chiral Lagrangean as done in 
Refs. [ill [H, [ItJ as a convenient starting point. The 
7r — a Lagrangean reads 

V(a)-^{U + Ui), (1) 

where U{x) — e lT ' 7T ^ 7r is the non-linearly transforming 
pion field and (,) represents the trace in isospin space. 
The scalar field is invariant under chiral transforma- 
tions 3 and the potential is chosen to have a minimum 
at o~ = f n . The sigma mass is then rr? a = V" {a)\ a= f^ , 
so that the physical scalar field is defined by the fluctu- 
ation around the vacuum expectation value, a = f x + s. 
}tc = 92.6 MeV denotes the pion weak-decay constant, 
ensuring the proper normalization condition of the pseu- 
doscalar fields. The vector mesons kinetic Lagrangeans 
are represented by Proca fields 

(2) 

and the kinetic nuclcon Lagrangian is 

The chirally invariant form of the meson-nucleon La- 
grangean can be looked up in Ref. [l5[. From the vac- 
uum expectation value of the scalar meson we get the 
nucleon mass Mn = g^NNf^ and the relevant nucleon- 
meson interaction vertices can be obtained from a chiral 



3 This is unlike the standard assignment of the linear sigma-model 
where one takes (u, 7?) as chiral partners in the (1/2, 1/2) repre- 
sentation of the chiral SU(2)n ® SU(2)l group. 
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Lagrangean [l_5j, |lg, [l7| and read 
ZA N 

C P NN = -9 p nnNt ■ p M 7 M iV 



JpNN 



2A 



N 



ZA N 



(4) 



Here and An is a mass scale which we take as An = 
3M/v/iV c with N c the number of colours in QCD. An 
overview of estimates of couplings from several sources is 
presented in appendix [XJ In the large N c limit the La- 
grangean simplifies tremendously since one has the fol- 
lowing scaling relations [4£| 4 



M N 
A N 

9ttNN 
fuNN 



N c , 
N° 

ly c > 

9aNN 

9pNN 



g^NN ~ fpNN 



'Nr. 



N' 



1/N C 



(5) 



The vector /tensor coupling dominance for u>/p is well 
fulfilled phcnomenologically (see appendix |A"| . Thus, in 
the large N c limit it is convenient to pass to the heavy 
baryon formulation by the transformation 



N(x) 



AMjyV-X 



B(x). 



(6) 



where B(x) is the heavy iso-doublet baryon field and v* 1 a 
four- vector fulfilling v 2 — 1, eliminating the heavy mass 
term [6(| [6lj]. Choosing v 1 * — (1,0) the meson- nucleon 
Lagrangean becomes 

£ = -g a NNsB^B + g uNN uj B^B 

+ g^NNB^aiTaBd'r + igZ-j* k Bt ai T a B& P ka . 

ZA N 

(7) 

In the large iV c -limit the contracted SU (4) algebra with 
the generators given by the total spin Si — J2A a f/^^ 
the total isospin 

T a = T,A T a/ 2 and the Gamow-Teller 
Xia = ^2,A a t T a I ^ operators is satisfied [5(| HJ). One 
could, of course, have started directly from the heavy- 
baryon Lagrangean, Eq. l[Tj). but the connection with chi- 
ral symmetry, in particular the relativistic mass relation 
Mn = ^o-jvat/tt would be lost. 



B. OBE potentials at leading N c 

From the heavy-baryon Lagrangean, Eq. ((TJ) , the calcu- 
lation of the NN potential in momentum space is straight- 
forward [l(| EH- However, passing to coordinate space 
is somewhat tricky since distributional contributions pro- 
portional to 6(x) and derivatives may appear. We discard 
them but just assuming that r > r c where r c is a short 
distance radial cut-off 5 . According to their increasing 
mass the leading N c contributions to the OBE potentials 
read 



<J\ ■ (72 - 



s 



V a {r) 
V P (r) 



12" 



„2 -m„r 
9a NN e 

47T r 

1 _ _ foNN ' 

Tl • To 

12 



4vr A N 



2(7i • rj 2 - 



■)2 



S 



12" 



9uiNN 



(m p r) 2 



(8) 
(9) 

(10) 
(11) 



4tt r 

where the tensor operator St2 = 3eri • xa% ■ x — o\ - 02 has 
been defined. Thus, the structure of the leading large 
A^-OBE potential has the general structure [52] 

V(r) = V c {r) +tx-T2 [01 • (7 2 W s {r) + S 12 W T (r)} . (12) 

Thus, we have as the only non-vanishing components 



V c (r) 
W s {r) 
W T (r) 



9oNN e 



9u,NN e 



4tt 

2 



1 9l 



NN 



47r r 

J- J o 



12 4tt A 2 n 

1 91nn m l e ' 
12 4tt A 2 n 7 

1 fpNN m 2 p e~ r 
12 4tt A^ 7 



r 



+ 6 " 4tt A2, r 



3 

m„r 



(m T r) 2 
3 



(m p r) 2 



At short distances we have 

9uNN ~ 9aNN 1 



Vfc(r) 
W s (r) 
Wr(r) 



47r r 

1 g^JVAf"^ + 2fpNN m2 p 1 

12 47TA2, r 

1 9rrN N — /pAW 1 

4 47TAI, r 3 ' 



(13) 

(14) 
(15) 

(16) 



4 There should be no confussion in forthcoming sections when we 
take N c = 3 and Ajv = Mn and the book-keeping becomes less 
evident. 



5 As discussed at length in Ref. [48l . l62l | these terms are ef- 
fectively inessential under renormalization of the correspond- 
ing Schrodinger equation via the coordinate boundary condition 
method, which will be explained shortly. 
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As we see, the potential is singular at short distances 
except for the very special value J p nn — 9irNN (see Ap- 
pendix [EJ . While the central Vc and spin Wg contri- 
butions present a mild Coulomb singularity, the tensor 
force component Wt is a more serious type of singular- 
ity, a situation appeared already for the simpler OPE 
potential fli ]. 



For the uncoupled spin-triplet channel we have 
~2R J(J + 1) 



v'Kr) 



vj{r) = p 2 vj{r) . (21) 



At very short distances we may neglect the centrifugal 
barrier and the energy yielding 



The OBE potential and ambiguities in the 
S-matrix 



We will show next that the S-matrix associated to the 
OBE potential is necessarily ambiguous, precisely be- 
cause of the short distance 1/r 3 singularity in the non- 
exceptional situation g-^NN ^ fpNN ■ The exceptional 
case, g V NN = fpNN will be treated in Appendix IBl We 
do so by proving that the standard regularity conditions 
for the wave function do not uniquely determine the solu- 
tion of the Schrodinger equation. Actually, at short dis- 
tances, i.e. much smaller than meson masses, r <C 1/jn,, 
the NN problem due to the OBE potential corresponds to 
the interaction of two spin-1/2 magnetic dipoles, namely 



V 2 * fe (f) + U dd (x)V k (x) = pH k {x) 



r <C 1 /m , 



(17) 



where the reduced dipole-dipole potential 6 is given by 

U dd (x) = MV dd (x) = 
R 

= (3fii • xa 2 ■ x — ax ■ 02) , (18) 



with R is a length scale and in our particular case 



±R = 



16ttA^ 



(9kNN 



fpNN) 



(19) 



the positive or negative sign depends on whether g^N n > 
fpNN or g-^NN < fpNN respectively. 

The above potentials become diagonal in the standard 
total spin 5 2 , parity Up, isospin T and total angular 
momentum J = L + S basis, so the states are labeled 
by the spectroscopic notation 2S+1 Lj. We remind that 
Fermi-Dirac statistics implies (— i) L + s '+ T = — 1. Thus, 
n ■ t 2 = 2T(T + 1) - 3 and a x ■ 02 = 25(5 + 1) - 3. 
For spin singlet states 5 = and J = L, the parity is 
natural Up = (— l)' 7 and one has S12 = 0. For uncoupled 
spin triplet states 5 — 1 one has J — L, natural parity 
Up = (— l) L and S12 = 2. For coupled spin triplet states 
5 = 1 one has L = J±l, unnatural parity Up — (— l) J+1 
and 



5l2 — 



2(J-1) 6^/7(7+1) 

2,7+1 2J+1 

2(7+2) 
2.7+1 



6^7(7+1) 
27+1 



(20) 



Note that we are not assuming here this potential at large dis- 
tances and so the standard long range problems of dipole scat- 
tering never appear. 



2 D 

-v'j(r)±— vj(r) =0, r<l/m,R,l/p. (22) 

The general solution can be written in terms of Bessel 
functions. Using their asymptotic expansions we may 
write at short distances 7 



v +,j( r ) 



v -,j( r ) 



r \ 3/4 r 
R 



/r\3/4 

\r) 



C\Ae~ 



C 2 pe 



C 2A e 4 



(24) 



Clearly in the repulsive case the regularity condition 
fixes the coefficient of the diverging exponential to zero, 
Cit? = 0, whereas in the attractive case both linearly in- 
dependent solutions are regular and the solution is not 
unique. In the case of the triplet coupled channel, we 
have for r <C 1/m, R, 1/p, i.e. neglecting centrifugal bar- 
rier and energy, the system of two coupled differential 
equations becomes 



-w»(r)) ± .3 



2(J-1) 6^7(7+1) 



2.7+1 



27+1 



/-' \ 6^/7(7+1) 2(7+2) 
27+1 27+1 




(25) 



This system can be diagonalized by going to the rotated 
basis 



vi,j(r) 
V2,j{r) 



2J+1 

7+1 
27+1 



7+1 
27+1 



J 



2,7+1 



w,j(r) 



-(26) 



7 The solutions of —y"(x) — y(x)/x s = are 

3 

x 1 



v/iYi(2/VI) 



003(77/4 + 2/^) + . . 



cos(7r/4 - 2/y-r) + . 



whereas the solutions of —y"(x) + y(x)/x' s = are 

v / 5Ji(2/V5) = — i —xie 2 /^+... 
2 v / 7r 



(23) 
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FIG. 1: The potentials Vc(r), Ws(r) and Wt{t) in MeV as a function of the distance (in fm). We include the effect of both 
exponential, Eq. (|31|l . and monopole Eq. (130[) form factors for A^jviv = 1300MeV and A^nn = 2000MeV. All other cut-offs 
are kept to A CT iviv = Ap^viv = A^jviv = 2000MeV. 



where the new functions satisfy 



AR 

v "A r ) T -Za v i,j( r ) 
-v'l.j{r) ± ^v 2 j(r) 



(27) 
(28) 



Note that here the signs are alternate, i.e. when one of 
the short-distance eigenpotentials is attractive the other 
one is repulsive and viceversa, and hence the type of solu- 
tions in Eq. (|24|) can be applied. This means that in gen- 
eral there will be solutions which are not necessarily fixed 
by the regularity condition at the origin, and thus the 
OBE potential does not predict the »!?— matrix uniquely. 
Instead, a complete parametric family of S— matrices will 
be generated depending on the particular choice of lin- 
early independent solutions, which are not dictated by 
the OBE potential itself. 

Thus, some additional information should be given. 
The traditional way is to introduce form factors to kill 
the singularity so that the regularity condition fixes the 
solution uniquely as we discuss in Section HTT1 Another 
way, which we discuss in the rest of the paper, is to fix 
directly the integration constants from data with or with- 
out form factors. As we will show this way of proceeding 
does not make much difference showing a marginal influ- 
ence of form factors (see Sect. I VII] ). 



Note that this assumes 1) Off-shell independence and 2) 
The form factor is accurately known. Standard choices 
are to take form factors of the monopole [l0| and expo- 
nential p} parameterizations 



1 mNN 



(<7 2 ) 



A 2 



i i 
q — m 



A 2 



(30) 
(31) 



T ™nn(<1 2 ) = ex P 
fulfilling the normalization condition r m ;yjv 

(m 2 ) = 1. 

These forms are so constructed as to have the same slope 
at small values of q 2 in the large cut-off expansion 

2 2 

T mNN (q 2 ) = 1 + A2 + 0(A- 4 ) . (32) 

so that the meaning for the cut-off is similar. In coor- 
dinate space this can be easily implemented for Yukawa 
potentials using 



yielding 



d 3 q 



(27t) 3 q 2 + m< 



-Ar 



[r m NN(q 2 )]' 



Y?° n (r) 



Airr 



Airr 



2A 



(33) 



(34) 



III. THE STANDARD APPROACH TO OBE 
POTENTIALS WITH FORM FACTORS 

A. Features of Vertex Functions 

A way out to avoid the singularities is to implement 
vertex functions in the OBE potentials corresponding to 
the replacement (q 2 = q$ — q 2 is the 4-momentum) 



NN 



(q) -> v mNN (q) [r mNN (q 2 )]' 



(29) 



which at short distances becomes finite, 



M-i^ + OCO, (35) 



which diverges linearly for A — > oo. The exponentially 
regularized Yukawa potential reads 



8irr 



8nr 



-Erf 



A 2 r — 4m 
2\/2A 



^ Erfc t^rJ ' (36) 
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where Erf and Erfc are the error function and comple- 
mentary error function respectively 8 . For Ar < 1 wc 
have the finite result 



-2m 2 /A 2 A m 

Erfc 

47T 



27T47T 



A 



0(r 2 ),(37) 



which diverges linearly for A — > oo. In the limit Ar ^> 1 
behaves as 



Yr(r) = 



Airr 



\A 2 r 2 e -2m 2 /A 2 

litAixr 2 



(38) 



\ NoFF(g* fflNN = 


19.7) — 


; \\ NoEF(g* fflNN = 


8.09) - 


\\\ MonFF (g* ffl NN = 


19.9) 


,\ M ° nFF (g*o)NN = 


3.05) 


\ E "P FF (g*oiNN = 


19.2) 


ExpFF (g*^ = 


2.87) .... 






>f \ 


= 1.3GeV - 


- I A a = A o. = A p 


= 2.0GeV - 



0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 
r[fm] 



and the distortion of the original Yukawa potential is 
much more suppressed in the exponential than in the 
case of monopole form factor. 

In any case we note the amazing feature that the form 
factors have a radically different effect on different com- 
ponents of the potential. While Vc and Ws with a mild 
~ 1/r short distance behaviour become finite, the tensor 
force behaving as Wt ~ 1/r 3 vanishes at the origin after 
due to the form factors, W£ wn (0) = W£ xp (0) = 0. This 
can be seen from the expression 



lim 



r^O J (27r) 3 



d 3 q ax ■qoi-q r . 2V] 



(2tt) 3 



q 2 + m 2 

[T m NN(q 2 )Y -m 2 Y x (0) 



d 3 q ! r , , - 



(39) 



which corresponds to take an angular average at short 
distances. This feature suggests that the impact of the 
tensor force at short distances should be small and looks 
clearly against the result of the short distance analysis 
outlined in Section [TT] where there is a strong mixing at 
short distances. As we will show in Sect. I VII within the 
renormalization approach there is no contradiction; phys- 
ical observables will naturally display a small mixing 9 . 

We show in Fig. Q] the potentials Vc(r), Ws(r) and 
Wt(t) in MeV as a function of the distance (in fm). We 
also include the effect of both exponential, Eq. (f3"Tj) . and 
monopole Eq. (f3"0"| form factors for A^nn = 1.3GeV and 
Attnn — 2GeV. All other cut-offs are kept to the val- 
ues A a NN = ApArjv = A^nn — 2GeV. As we see, the 
distortion of the tensor component due to the strong 
form factor takes place already at r ~ lfm for soft- 
est cut-off A^nn = 1.3GeV. The key issue here is to 



8 They are defined as 

Erf(je) = 1 - Erfc(z) = — f 
VV Jo 



dt e 



[l + Oiz- 1 )} 



9 In other words, the counterterm structure is not of the naive 
form suggested by Eq. I|39| l but a more general one including the 
tensor operator 5i2 . See also the discussion in Sect. IIVI 



FIG. 2: The 1 So potential Vl s (r) in MeV as a function of 
the distance (in fm) for the different scenarios with large and 
small bj— couplings. We include the effect of both exponential, 
Eq. (|31[) , and monopole Eq. (|30[) form factors for A^-jvat = 
1300MeV. All other cut-offs are kept to A^nn = A p nn = 
A^jviv = 2000MeV. 



decide whether this distortion represents a true physi- 
cal effect, rather than a mere artifact of the regulariza- 
tion. This boils down to determine if one can visualize 
finite nucleon size effects when the probing wavelength 
is not shorter than 0.5fm < r < lfm. The fact that the 
monopole and exponential parameterizations agree down 
to r ~ 0.5fm but differ from the bare unregularized po- 
tential suggests that one could look for a true physical 
effect based on model independent distortions in the re- 
gion slightly above 0.5fm. This point will be analyzed 
further in Section IVlII 

Finally, note that the multiplicative manner in which 
form factors are introduced, although it looks quite nat- 
ural, does build in correlations which may not reflect the 
real freedom one has in general; a dynamical calculation 
need not comply to this factorization scheme. 



B. The problem of short distance sensitivity vs 
spurious bound states 

The advantage of using vertex functions is that they 
make the OBE non-singular at short distances. As a con- 
sequence, the choice of the regular solution determines 
the solution uniquely. In this section we analyze critically 
the use of form factors which are customarily employed 
in NN calculations based on the OBE potential. We will 
see that for natural choices of meson-nucleon parameters 
(see Appendix [A"]) , the NN potential displays short dis- 
tance insensitivity and at the same time spurious deeply 
bound states. However, if we insist on not having spuri- 
ous bound states the resulting description is highly short 
distance sensitive. 

As we have mentioned, NN scattering in the elas- 
tic region below pion production threshold involves CM 
momenta p < p ma x — 400 MeV. Given the fact that 
\/m u r~j 1/rrip ~ 0.25fm <C l/pmax — 0.5fm we ex- 
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pect heavier mesons to be irrelevant, and uj and p to be 
marginally important, even in s- waves, which are most 
sensitive to short distances. This desirable property has 
not been fulfilled in the traditional approach to OBE 
forces. In order to illustrate this, we consider the 1 S$ 
channel, where the potential (without form factor) is 



Vt So (r) = Vc(r)-3W s (r) 



2 2 - 

9irNN m -K e 

2 p -m u r 
9uNN e 

47T r 



pNN 



4tt 

ml e - 



SttM 2 n 



(40) 




NoFF(g* fflNN = 19.7) 
NoEF (g* mNN = 8.09) 

MonFF(g* mNN = 19.9) 

MonFF(g* mNN = 3.05) —. 

ExpFF(g* (1)NN = 19.2) 

ExpFF(g* (oNN = 2.87) .... 



0.5 



I 

r[fm] 



1.5 



We take 



Tfi-x = 138MeV, M N = 939MeV, m p = 770MeV, 
m u) = 783MeV and g^NN = 13.1, which seem firmly 
established, and treat m a , g a NN and g^NN and f P NN 
as fitting parameters. To see the role of vector mesons 
we note the redundant combination of coupling constants 
9wNN- fpNN m l/ (2 m n) which appears in the 1 S poten- 
tial when we take m p = m w , a tolerable approximation 
within the present context. To avoid unnecessary strong 
correlations we define the effective coupling 



FIG. 3: Zero energy wave function for the singlet pn 1 So 
channel as a function of distance (in fm) and for the different 
scenarios with large and small lu— couplings. We include the 
effect of both exponential, Eq. (|31[) , and monopole Eq. (|30[) 
form factors for A-^nn = 1300MeV. All other cut-offs are kept 
to A a MN = Apjvjv = Aojjvjv = 2000MeV This wave function 
goes asymptotically to uo(r) — > 1 — r/ao with ao — — 23.74fm 
the scattering length in this channel. The zero at about r = 
0.5fm signals the existence of a spurious bound state. 



9uNN 



f, 



pNN m "p 



2M% 



(41) 



Natural values for the coupling constants g^NN = 9—10.5 
a.ndf pNN = 14-18 imply g* NN = 0-7. In Fig. [2]the po- 
tential without and with monopole and exponential ver- 
tex functions is depicted for several values of g t 



As 

we see, the differences start below lfm where the stan- 
dard short distance repulsive core is achieved by large 
and unnatural values of ff^Af' an< ^ n °t so mucn depend- 
ing on the form factors. On the other hand, if we use the 
regularized potential at r = and take the natural 
values for the coupling constants g^NN — 9 — 10.5 and 
fpNN = 14—18 the potential at the origin becomes 



Vi So (0) = -(1000- 3000)MeV. 



(42) 



which is huge and attractive. The number of states is 
approximately given by the WKB estimate 



T JO 



y/-MVi So (r)dr. 



(43) 



which yields numbers around unity. In fact the potential 
accommodates a deeply bound state, at about 



In the standard approach the scattering phase-shift 
5o(p) is computed by solving the (s-wave) Schrodinger 
equation in r-space 

-u';(r)+M N V(r)u p (r) = p 2 u p (r) (45) 

sm(pr + 6g(p)) 

Up{r) -► , (46 

sm d {p) 

with a regular boundary condition at the origin 

«p(0) = (47) 

This boundary condition obviously implies a knowledge 
of the potential in the whole interaction region, and it 
is equivalent to solve the Lippmann-Schwinger equation 
in p-space. In the usual approach [ToL [l3| everything is 
obtained from the potential assumed to be valid for < 
r < oo. In practice, and as mentioned above, strong form 
factors are included mimicking the finite nucleon size and 
reducing the short distance repulsion of the potential, 
but the regular boundary condition is always kept. One 
should note, however, that due to the unnaturally large 
NN 1 So scattering length (ao ~ — 23fm), any change in 
the potential V — > V + AV has a dramatic effect on ao, 
since one obtains 



E F 



-(500- 2000)MeV 



(44) 



This state does not exist in nature and should clearly be 
ruled out from the description on a fundamental level. 
On the other hand, we do not expect such state to influ- 
ence the low energy properties below the inelastic pion 
production threshold -Ecm — 175MeV in any significant 
manner. 



Aq = ttpMjv / AV(r)u (rydr 



(48) 



a quadratic effect in the large ao- This implies that po- 
tential parameters must be fine tuned, and in particular 
the short distance physics. To illustrate this we make 
a fit the np data of Ref. [HI]. The results using the 
OBE potential without or with strong exponential and 
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TABLE I: Fits to the 1 So phase shift of the Nijmegen group [63l ] using the OBE potential without or with strong exponential 
and monopole form factor. We take m = 138.03MeV, and g^NN = 13.1083 [sH and m p = = 770MeV and fit rria, g a NN 
and g^NN- We use A^jvat = 1300MeV and A^nn = A p jvjv = A^jvat = 2000MeV. Eb represents the energy of the (spurious) 
bound state when it does exist. 





r c (fm) 


m CT (MeV) 


gaNN 


9uNN 


X 2 /DOF 


ao (fm) 


ro (fm) 


Ss(MeV) 


r(g 2 ) = i 





547.55(4) 


13.559(8) 


19.68(2) 


0.869 


-23.742 


2.702 




r(? 2 ) = i 


0.1 


500.9 (5) 


9.61(1) 


8.09(2) 


0.484 


-23.742 


2.504 


-638 







552.57(5) 


13.78(2) 


19.21(4) 


0.664 


-23.741 


2.703 




r(9 2 )=r axp (?2) 





525.1(1) 


10.41(1) 


2.9(1) 


0.213 


-23.740 


2.698 


-578 


r( q 2 ) = r mon ( g 2 ) 





551.7(1) 


13.99(1) 


19.978 (11) 


0.971 


-23.741 


2.707 




Y{q 2 ) = r mon ( 9 2 ) 





532.5(2) 


10.81(1) 


3.04(3) 


0.241 


-23.739 


2.696 


-597 



monopole form factor 10 are presented in Table Q] In all 
cases we have at least two possible but mutually incom- 
patible scenarios. An extreme situation corresponds to 
the case with no form factors 11 . The small errors should 
be noted, in harmony with the fine tuning displayed by 
Eq. (|48p and the corresponding couplings and scalar mass 
are determined to high accuracy but turn out to be in- 
compatible. This is just opposite to our expectations 
and we may regard these fits, despite their success in de- 
scribing the data, as unnatural. The ambiguity in these 
solutions are typical of the inverse scattering problem, 
and has to do with the number of bound states allowed 
by the potential. Actually, this can be seen from Fig. [3] 
where the zero energy wave function is represented. Ac- 
cording to the oscillation theorem, the number of interior 
nodes determines the number of bound states. Thus, the 
larger values of S^aw correspond to a situation with no- 
bound states since uo(r) does not vanish, whereas for the 
smaller g^, NN values one has a bound state as Uo(r) has 
a zero, which energy can be looked up in Table Q] Of 
course, such a bound state does not exist in nature and 
it is thus spurious. On the other hand they always take 
place at more than twice the maximum energy probed 
in NN scattering, Eqm = 175MeV, and we should not 
expect any big effect from such an state. Note that de- 
spite the net repulsive w-vector and attractive p-tensor 
couplings, the total potential would not be repulsive at 
short distances with or without form factors in the solu- 
tion with natural couplings and a spurious bound state; 
the net short distance repulsion comes about only in the 
solution with unnaturally large coupling (see Fig. [5]). 

From the table U one can clearly understand the usu- 



ally too large values of the g^NN coupling constant as 
compared to those from SU(3) symmetry g W MN ~ 9 or 
from the radiative decay u> — > e + e _ yielding g^NN = 
10.2(4). Using the definition of g* NN , Eq. !j4"I|). we get 
for fpNN = 14—18 large values of g^NN = 20 — 22 for the 
case with no bound state, whereas more natural values 
g^NN = 8.5 — 10.5 are obtained for the case with one 
(spurious) bound state. 



IV. BOUNDARY CONDITION 
RENORMALIZATION AND ULTRAVIOLET 
COMPLETENESS 



According to the discussion of Sect. Ill Cl the short dis- 
tance 1/r 3 singularity of the OBE potential makes the 
solution ambiguous, and thus there is a flagrant need of 
additional information not encoded in the potential it- 
self. Of course, once we realize the freedom of choosing 
suitable linear combinations of independent solutions, we 
may question how general this choice can be, even if the 
potential is not singular. In this Section we derive con- 
straints on the short distance boundary condition. As 
mentioned already, we work with energy independent po- 
tentials. In this section we show what this requirement 
implies for the renormalization program. Using the po- 
tential of Eq. (TT2")) we solve the Schrodinger equation, 



A I 



V 2 * fe (f) + V{x)V k (x) = E k V k (x) , (49) 



where ty(x) is a spin-isospin vector with 4x4=16 compo- 
nents, which usually satisfies the out-going wave bound- 
ary condition, 



10 In this particular channel the regularity condition, Eq. 1471 deter- 
mines the solution completely since the potential without vertex 
functions V15 (r) ~ 1/r is not singular at short distances in the 
sense that lim r ^ 2fj,\V(r)\r 2 = 00 [2Cll2ll. 

11 We find strong non-linear and well determined correlations have 
been found making a standard error analysis inapplicable. In 
this situation we prefer to quote errors by varying independently 
the fitting variables g a NN, W CT and g^ NN until A% 2 = 3.53 as 
it corresponds to 3 degrees of freedom. 



*fc(x) 



ik-x 



p ikr 

f(k',k)— 

r 



s,m s 
/Ct,mt 



(50) 



with f(k , k) the quantum mechanical scattering matrix 
amplitude and Xt'mt a 4x4 total spin-isospin state. We 
apply a radial cut-off r c and consider that the local po- 
tential V(x) is valid for the long distance region r > r c . 
The precise form of the interaction for the short distance 
region r < r c is not necessary as the limit r c — > will 
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be taken at the end. To fix ideas we assume an energy 
independent non-local potential, as we expect genuine 
energy dependence to show up as sub-threshold inelastic 
(e.g. pion production) effects. Any distributional terms 
~ 5(x) arising from the long distance potential V{x) are 
necessarily included in the inner region, r < r c . The 
inner wave function $> k (x) satisfies 

- -^V 2 $ fc (x) + J d 3 x'V(x',x)$ k (x') = E k $ k (x) ,(51) 

and will be assumed to vanish at the origin. Using stan- 
dard manipulations and the Green identity we get for the 
inner and outer regions 

(E p — E k ) I d?x$l{x)%{x) 

J r<r c 

v$i(f)$ p (f) - $[. (f)v$,(f) 



dS 



,(52) 



and 



[E p — E k ) f d 3 x^l(x)^ p (x) 

J r~>r c 



= - dS 



- ^(x)W p (x) 



(53) 



respectively where the difference in sign form the inner to 
the outer integration comes from opposite orientations in 
the integration surface. Clearly, orthogonality of states 
in the whole space for different energies, 



d\$>i{x)%{x) 



d z x^\(x)^p{x) =0,(54) 



can be achieved by setting the general and common 
boundary condition, 



d r $ p (xr c ) = Lp(xr c )$ p (xr c ) 
d r ^ p (xr c ) = Lp(xr c )^ p (xr c ) . 



(55) 



Here, L p (xr c ) is a self-adjoint matrix which may de- 
pend on energy, and may be chosen to commute with 
the symmetries of the potential V(x) 12 . Deriving with 
respect to the energy the inner boundary condition, 
Eq. (52]), i.e. taking E p = E k + AE and ® p (x) = 
$ fc (x) + AEd$ p (x)/dE , we get 



dx^p{xr c ) d -^^%(xr c ) 



r 2 dr I dx$Uxr)$ p (xr) 



Mr 2 c Jo 



(56) 



12 In practice this would mean taking 

L(xr c ) = L c {r c ) + ri ■ r 2 [Is(r c )ffi ■ cr 2 + L T (r c )S 12 ] 
which implies at most only three counterterms for all partial 



where we see that MdL p (xr c )/dE 



The important 



issue here is that regardless on the representation at short 
distances, the boundary condition must become energy 
independent when r c — > 0, namely 



]im 9L p (xr c ) = 

r c ^0 dE 



provided one has 



lim 



dr I dx$p(xr)$ p (xr) = . 



(57) 



(58) 



Thus we may take a fixed energy, e.g. zero energy, as a 
reference state. 



lim 



L p (xr c ) 



lim Lo(xr c ) . 



(59) 



The condition of Eq. (|58[) is the quite natural quantum 
mechanical requirement that the contribution to the to- 
tal probability in the (generally unknown) short distance 
region is small. This is the physical basis of the renormal- 
ization program which corresponds to the mathematical 
implementation of short distance insensitivity and which 
we carry out below for the OBE potential. It should 
be noted that this requirement depends on the poten- 
tial. The condition of Eq. (f58"|) implies that in the limit 
r c — ► one must always choose a normalizable outer so- 
lution ^ k (x) at the origin and the boundary condition 
must be chosen independent on energy. Note that en- 
ergy dependence would be allowed if the cut-off was kept 
finite, and still the requirement of orthogonality in the 
whole space could be fulfilled for an interaction charac- 
terized by a non-local and energy independent potential 
in the inner region. This simultaneous disregard of both 
non-local and energy dependent effects was advocated 
long ago by Partovi and Lomon Q on physical grounds, 
and as we see, it is a natural consequence within the 
renormalization approach. 

The renormalization procedure is then conceptually 
simple since any energy state with given quantum num- 
bers can be chosen as a reference state to determine the 
rest of the bound state spectrum and scattering states. 
For, instance using a bound state (the deuteron) ^(x), 
at long distances (see Sect. IVII for details) 



A, 



1 + Ts s " 



y sm s 



(60) 



we integrate in the deuteron equation 



1 



-^-V 2 ^(f) + V($)* d {x) = -^ d (x) , (61) 



M 



M 



and determine the short distance boundary condition ma- 
trix L[xr c ) from 



d r ^d(xr c ) = L(xr c )^ d (xr c ) 



(62) 



Then, using the same boundary condition for the finite 
energy state, 



d r ^ k (xr c ) = L(xr c )^ k (xr c ) , 



(63) 
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we integrate out the finite energy equation (|49|) whence 
the scattering amplitude may be obtained. In this man- 
ner the deuteron binding energy defines the appropriate 
self-adjoint extension spanning the relevant Hilbert space 
in the 3 Si — 3 D\ channel. Rcnormalization is achieved 
by taking the limit r c — ► at the end of the calculation. 
The conditions under which such a procedure is mean- 
ingful will be discussed below for the particular partial 
waves under study, but a fairly general discussion can be 
found in Refs. 0, [46|, [62j . Relevant cases for chiral po- 
tentials where this condition turned out not to be true 
are discussed in Ref. (62J. We will also encounter below 
a similar situation in our description of the deuteron and 
the 3 Si — 3 D\ channel. As a consequence of the previous 
limit the completeness relation reads 



energy plane. This restriction complies to causality, and 
implies in particular the fulfillment of Wigncr inequalities 
as was discussed for the 1 Sq channel in Ref. [62j 14 . 



V. THE SINGLET CHANNEL 

A. Equations and boundary conditions 

The wave function in the pn CM system can be 
written as 



/,/3 h. 
— ^ k {^{(x') + J2 *n(J)*J,( 
! ~" ' £„<0 



') = S(x-x')l. 



Besides the deuteron, the sum over negative energy states 
contains most frequently spurious bound states, and for 
the singular potential such as the one we are treating 
here there are infinitely many. They show up as oscilla- 
tions in the wave function at short distances, and are a 
consequence of extrapolating the long distance potential 
to short distances. On the other hand, from the above 
decomposition one may write a dispersion relation for the 
scattering amplitude 13 of the form 



f{k',k) 



4tt 



E — E n 



E n <0 

M f d 3 q (%\V\* q ){* q \V\k) 



4-7T 



(2tt) 



E-q 2 /M 



(65) 



where /s(fc',fc) is the Born amplitude and the physical 
and spurious bound states occur as poles in the scat- 
tering matrix at negative energies E = E n and the dis- 
continuity cut along the real and positive axis is given 
by the second term only. Clearly, the influence of these 
spurious bound states poled is suppressed if their energy 
E n <C E c i- Given the fact that these states do occur in 
practice it is mandatory to check their precise location 
to make sure that they do not influence significantly the 
calculations, or else one should study the dependence of 
the observables on the short distance cut-off r c starting 
from a situation where it is small but still large enough 
as to prevent the occurrence of spurious bound states. It 
should be noted, however, that in no case can the spuri- 
ous states occur in the first Riemann sheet of the complex 



1 



Air r 



(66) 



with the total spin s = and m s = 0. The function u(r) 
is the reduced S- wave function, satisfying 



(64) 



u'p(r) + Ui So (r)u p (r) = p 2 u p (r) 



where one has 



Ui So = M(V C - 3W S ) 



(67) 



(68) 



At short distances the OBE potential behaves as a 
Coulomb type interaction, 



±- 



1 

Rr 



(69) 



where 



1 _ M 
R ~ 4^ 



2 

9u>NN 



f2 2 
JpNN m p 

2M 2 N 



9<tnn Unn 70) 



Here, f^NN = g^N ntti-k / {^M n) . The repulsive or attrac- 
tive character of the interaction depends on a balance 
among coupling constants. The short distance solution 
can be written as a linear combination of the regular and 
irregular solution at the origin 



u p (r) -> ci(p) + c 2 (p)r/R 



(71) 



where in principle the arbitrary constants ci (p) and C2 (p) 
depend on energy To fix the undetermined constants we 



13 This is done by using the Lippmann-Schwinger equation in the 
form T = V + VGV with G = (E - H)' 1 , and normaliza- 



tion (k\x) 
-M/(4ir)(k'\T(E)\k). 



and { "J/ | a?) = ^^(x) whence f(k',k) 



14 Causality violations, i.e. poles in the first Riemann sheet of 
the complex energy plane are easy to encounter (see e.g. [621), 
particularly with energy dependent boundary conditions. A 
prominent example is an s-wave without potential and having 
u' p (0)/u p (0) = pcot Sq(p) = — 1/qo + r 0P 2 /2+i | 2P 4 , which for the 
1 So-channel values of parameters ao = — 23.74fm, r*o = 2.75fm 
and V2 = — 0.48fm 3 yields besides the well-known virtual state 
in the second Riemann sheet E v = — 0.066MeV a spurious 
bound state at Eg = — 18.37MeV and an unphysical pole at 
E = 128.88 ± i46.45MeV. However, finite cut-offs and energy 
independent boundary conditions are guaranteed not to exhibit 
these problems, while some spurious bound states may be re- 
moved. 
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FIG. 4: Left: A\ 2 = 2.3 confidence level ellipse (corresponding to 68% for 2 variables) in the g a NN — m a plane without vector 
mesons q^nn = f P NN = 0. Right: Renormalized OBE 1 So pn phase shifts (in degrees) as a function of CM momentum. Data 
from 16311. 




FIG. 5: X 2 / D OF (left panel), gaNN (middle panel) and m a (right panel) as a function of the effective coupling constant 
9Znn = \JoIjnn fpNN m 'p/^^N wnen a fit of the leading N c contributions to the OBE (ct + ty + p + to) potential is considered. 



impose orthogonality for r > r c between the zero energy 
state and the state with momentum p and get 



u' p (r c )u k (r c ) - u' k (r c )u p (r c ) 

= (k 2 -p 2 ) / u k (r)u p {r)dr = 0,(72) 



B. Phase shifts 

For a finite energy scattering state we solve for the 
OBE potential with the normalization 



u p {r) 



sin(pr + 6 (p)) 
sin S (p) 



(74) 



Taking the limit r c — > implies the following energy 
independent combination [45| 



ci(p) = CiW = ci(0) 
ca(p) c 2 (fc) c 2 (0) 



(73) 



leaving one fixed ratio which can be determined from e.g. 
the zero energy state or any other reference state. 



with So(p) the phase shift. For a potential falling off 
exponentially ~ e - m ^ r a t large distances, one has the 
effective range expansion at low energies, \p\ < m w /2, 



pcotd {p) = -- + \r p 2 + V2P 4 " + . . . (75) 
ao 2 



with ao the scattering length and r$ the effective range. 
The phase shift is determined from Eq. (|74p . Thus, for 
the zero energy state we solve 



<(r) + Uk So (r)«o(r) = 0. 



(76) 
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with the asymptotic normalization at large distances, ob- 
tained from Eq. (|7"4"1) . 



Uo(r) 



1 



r 



(77) 



In this equation ao is an input, so one integrates in 
Eq. (|76|) from infinity to the origin. Then, the effective 
range defined as 



r = 2 



dr 



1-— -Mr) 2 
a , 



(78) 



can be computed. From the superposition principle of 
boundary conditions 



u (r) = uo, c (r) — u 0iS (r) 

a 



(79) 



where uo,cM — > 1 and uo, s (r) — > r correspond to cases 
where the scattering length is either infinity or zero re- 
spectively. Using this decomposition one gets 



r 



, B c 

a «5 



(80) 



where 



c 

Q 
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m o( r c) 
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r c [fm] 



.4 
B 

C 



2 / dr(l - < c ) , 
Jo 

poo 

-4 / dr(r - u 0:C uo,s) , 
Jo 

/>oo 

2 / dr(r 2 - < s ) , 
Jo 



(81) 
(82) 
(83) 



depend on the potential parameters only. The interesting 
thing is that all dependence on the scattering length ao is 
displayed explicitly by Eq. §5U\ ) . To determine the phase 
shift 5o(p) one proceeds as follows. From Eq. (1771) an d 
integrating in Eq. (|76p one determines ci(0) and C2(0) 
and uses Eq. ([73"]) to determine the ratio c\ (p) / (p) and 
integrates out Eq. (|67[) matching Eq. (|74|) . This way the 
phase shift 5o(p) is determined from the potential and the 
scattering len gth as independent parameters. As it was 
shown in Ref. [48| this procedure is completely equivalent 
to renormalize the Lippmann-Schwinger equation with 
one counterterm. 



C. Fixing of scalar parameters 

In this work we will fix our parameters in such a way 
that the 1 5'o phase shift is reproduced. This has the ad- 
vantage that the scalar meson parameters are determined 
for the rest of observables. Thus, fixing the scattering 
length ao = — 23.74fm and the OPE potential parame- 
ters g^NN = 13.1 and = 138.04MeV we fit g a NN and 
m a to the 1 Sq phase shift of the Nijmegen group [63j . 



FIG. 6: Short distance cut-off r c (in fm) dependence of the 
1 So phase shift fitting procedure. \ 2 /DOF (upper panel) and 
m a (lower panel). The bands represent the error estimate for 
two fixed cut-off values: a finite cut-off r c — 0.4fm and the 
renormalized case r c — Ofm. 



In the absence of vector meson contributions, i.e. taking 
9ujNN = fpNN = the fit yields 



g aNN = 9(1) m a = 501(25)MeV 



(84) 



with a x 2 /DOF = 0.13. As we see from Fig. 2J there 
is a large, in fact linear, correlation, between the scalar 
coupling and mass, while the fit is quite good as we can 
see. For comparison we also show the result with OPE 
which, despite reproducing the threshold behaviour does 
a poor job elsewhere. We quote also the effective range 
values from the universal low energy theorem, 



r = 1.3081 



4.5477 5.1926 



o 



id 



1.5089fm, 



v'o 



2.4567- 
2.6989fm 



5.5284 5.7398 



(tt) 



(tt + ct) 



(85) 



where the corresponding numerical values when the ex- 
perimental ao = — 23.74fm have also been added. 

It is interesting to analyze the dependence of the fit- 
ted scalar parameters on the short distance cut-off ra- 
dius, r c . A priori we should see the a exchange for 
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r c < l/m a = 0.4fm. From Fig. [6] we see the masses and 
the couplings providing an acceptable fit x 2 /DOF < l for 
which a reliable error analysis may be undertaken. As we 
see this happens for r c < 0.6fm and two stable plateau 
regions yielding two potentially conflicting central m a 
values. An error analysis both at a finite cut-off value 
r c = 0.4fm and the renormalized cut-off limit r c = Ofm. 
gives two overlapping and hence compatible bands. This 
shows that in this case the data do not discriminate below 
r c = 0.5fm. Much above that scale, the a meson becomes 
nearly irrelevant, as the coupling becomes rather small. 

Alternatively, we may treat the cut-off itself as a fitting 
parameter. To avoid the large m a — g a NN correlations 
displayed in Fig. [4] we fix the coupling constant to its 
central value g a NN = 9.1 and get then r c — O-lOj^Q Q^fm 
and m a = 500(3)MeV. This shows that removing the 
cut-off is not only a nice theoretical requirement, but 
also a preferred phenomenological choice. 

To analyze now the role of vector mesons we note, as al- 
ready discussed in Section UlI Bi the redundant combina- 
tion of coupling constants g^ NN — fpNN 171 ^ which 
appears in the 1 5'o potential when we take m p — m u . 
We thus define the effective coupling g^ NN defined in 
Eq. (|4ip. This combination is responsible for the re- 
pulsive contribution to the potential in the 1 Sq channel. 
From typical values of the couplings g^NN = 9 — 10.5 and 
fpNN = 15 — 17 we expect g^ NN to be effectively small. 
We show in Fig. ([5]) the corresponding x 2 /DOF as well 
as the readjusted scalar mass m a and coupling g a NN as 
a function of the effective combination of coupling con- 
stants, <?* NAT- As we see, the fit is rather insensitive but 
actually slightly worse than without vector mesons when 
their contribution is repulsive. Thus, we will fix this ef- 
fective coupling to zero which corresponds to take 



2 

9uNN 



f, 



pNN m 'p 



2M 2 N 



(86) 



This choice has the practical advantage of fixing g a NN 
and m a to the values provided in Eq. (|84p also when 
the leading N c vector meson contributions are included. 
Moreover, it is also phenomenologically satisfactory as we 
have discussed above. In Sect. IVII we will also see that 
deuteron or triplet 3 Si — 3 D\ do not fix the deviations 
from the relation given by Eq. 



D. Discussion 

The linear g a NN — "V correlation can be established 
solely by requiring that the effective range of the Ni- 
jmegen group ro = 2.67fm or any other be repro- 
duced [65(1 . Actually, Eq. (f8"4")l . yields the combination 
C a = g 2 aNN /ml = 331(50)GeV" 2 which is fixed by the 
effective range and not by the scattering length. This 
is in contrast with the resonance saturation viewpoint 
adopted in Ref. [66[ where this combination fixes the scat- 
tering length. 




t [fin] 



FIG. 7: Zero energy wave function for the singlet pn 1 5o 
channel as a function of distance (in fm) . This wave function 
goes asymptotically to Uo(r) — > 1 — r/ao with ao = — 23.74fm 
the scattering length in this channel. 



Furthermore our calculation shows that an accurate 
fit without explicit contribution of the vector mesons is 
possible. In particular, our potential does not exhibit any 
repulsive region. This is in apparent contradiction with 
the traditional view point that the w-meson is responsible 
for the short range repulsion of the nuclear force. 

To understand this issue we plot in Fig. [7] the zero en- 
ergy wave function obtained by integrating in with the 
physical scattering length ao- As we see, there appear 
two zeros indicating, according to the oscillation theo- 
rem, the existence of two negative energy spurious bound 
states. To compute such a state we solve Eq. (|57| with 
negative energy E B = —j B /M, for an exponentially de- 
caying wave function, ugfr) — * A B e~ lBT (normalized to 
one), and impose orthogonality to the zero energy state, 
namely 



u u B - u u B 



0. 



(87) 



from which 75 can be determined. A direct calculation 
yields E B1 = -777MeV and A B \ = 15.64fm _1/2 and 
E B2 = -11077MeV and A B2 = 27.43fnT 1/2 . If we re- 
gard the scattering amplitude as a function of energy in 
the complex plane, these spurious bound state energies 
are beyond well the maximum CM energy we want to 
describe in elastic NN scattering, Eqm < 175MeV, and 
so have no practical effect on the scattering region. The 
appearance of spurious bound states in EFT approaches 
are commonplace; one must check that they are beyond 
the considered energy range. 

In order to discuss this point further we may try several 
ways of removing the unwanted poles and to quantify the 
effect on the results. Unitarity implies the usual relation 
between the partial wave amplitude and the phase shift 



[fo(p)} =p cot 6 {p) -ip 



(88) 



Actually, the contribution of a negative energy state to 
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the s-wave scattering amplitude is a pole contribution 



VI. THE TRIPLET CHANNEL 



fo(p)\i 



A 2 



1 



A 2 



M E 



P 



(89) 



A simple way of subtracting such a bound state without 
spoiling unitarity and preserving the value of the ampli- 
tude at threshold /o(0) = F o (0) = — «o is to modify the 
real part of the inverse amplitude as follows, 



1 



1 



Fo(p) foip) 



P 
A 2 

B 



(90) 



which has no pole at E = —\Eb\, since Fo(ij E ) = 
A 2 B 1 7^ . Using the relation between amplitude and phase 
shift Fq(j>) = l/(pcot A (p) ~ ip) we get the modified 
phase shift, 



pcot A Q (p) = pcotS (p) + , (91) 
which corresponds to a change in the effective range 



Ar °' B = I| 



(92) 



For the values of the two spurious bound states we get 
Ar'ols = 0.008, 0.002fm, a tiny amount. The change in 
the phase shift never exceeds 0.1°. Of course, this is not 
the only procedure to remove spurious bound states, but 
the result indicates that the effect should be small. 

Another practical way to verify this issue is to study 
the influence of changing the cut-off r c from the lowest 
value not generating any spurious bound state and the 
origin, corresponding to look for u Q (a) = 0. This point is 
clearly identified as the outer zero of the wave function, 
which takes place at about a = 0.5fm. Thus, if we choose 
r c = a, there will not be any bound state. For this partic- 
ular point, the orthogonality of states, Eq. (|72|) . implies 
that u p {a) = 0, resembling the standard hard core pic- 
ture, if we assume u p (r) = for r < a. Thus, at this r c 
our method would correspond to infinite repulsion below 
that scale. In other words, the boundary condition does 
incorporate some effective repulsion which need not be 
necessarily visualized as a potential. The advantage of 
using a boundary condition is that we need not require 
modelling nor deep understanding on the inaccessible and 
unknown short distance physics. 

The contribution to the effective range from the ori- 
gin to the "hard core" radius a is ~ 0.04fm, while 
the change in the phase shift at the maximum energy 
due to the inner region < r < a is A<5o = 6° to be 
compared with the error estimate A<5n_ = 0.7° from the 
PWA analysis of the Nijmegen group [14j or the A^o = 2° 
from the corresponding high quality potentials [63| . If we 
identify this hard core radius to the breakdown scale of 
the potential, these differences might be interpreted as a 
systematic error of the renormalization approach for our 
OBE potential and, as we see, they turn out to be rather 
reasonable. 



A. Equations and boundary conditions 

The 3 Si — 3 Di wave function in the pn CM system can 
be written as 

1 



w(r) 



u(r)a p 



(3cr p 



cr„) 



X', 



(93) 



0, ±1 and o~ p and 



with the total spin s — 1 and m, 
o~ n the Pauli matrices for the proton and the neutron 
respectively. The functions u(r) and w(r) are the reduced 
S- and D-wave components of the relative wave function 
respectively. They satisfy the coupled set of equations in 
the 3 Si — 3 Di channel 



u"{r) + U3 Sl (r)u(r) + U El (r)w(r) 



-w"(r) + U El {r)u(r) + 



Us Dl (r) + 



w(r) 



MEu(r) , 
MEw(r) . 



(94) 



with Uss 1 (r), Ue 1 (t) and U^ Dl {r) the corresponding ma- 
trix elements of the coupled channel potential 

Us Sl = M(Vc-3W s ), 

U El = -6V2MW T , 

U3 Dl = M(Vc-3W s + 6Wt). (95) 
At short distances one has the leading singularity 
U, Sl = Oir- 1 ), 
4V2R 



U El 



+ 0(r- 1 ) 



where 



±R = 



a 2 ~ f 2 



(96) 



(97) 



32wM N 

This is very similar to the pure OPE case treated in 
Ref. but with the important technical difference that 
for fpNN < 9-kNN and / p nn > 9tiNN there is a turn-over 
of repulsive-attractive eigenchannels since the effective 
short distance scale R changes sign. Thus, we must dis- 
tinguish two different cases 15 . At short distances we 
have for g^NN > fpNN the plus sign in Eq. (|97|) yielding 



u A {r) 



u R {r) = --L u (r) + ^w{r), (98) 



15 The exceptional case, g^jviv = fpNN corresponds to a regular 
potential and will be treated in Appendix FBI 
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whereas for g n NN < fpNN the minus sign in Eq. (f97|) is 
taken and the solutions are interchanged 



UR(r) 



u A (r) = --L u (r) + y~u;(r), 

yielding an attractive singular potential XJ a 
for ua and Ur — > 8R/r 3 for ur, which solutions are 



(99) 



-4R/r 3 



ur{t) 



/r\3/4 

\r) 



(100) 



UA 



M - (5 



7?/ 



x 3/4 



C 



iAe 



C 2A e il 



The constants Ci/j, C 2 r, C\a and depend on both 
7 and 77 and the OBE potential parameters. As it was 
discussed in Ref. [44j we must define a common domain 
of wave functions to define a complete solution of the 
Hilbert space in this 3 <Si — 3 D\ channel. This is achieved 
taking 



u R (r) - C R {i) (£) 
«A(r) -> CU( 7 ) (£) 



. ,R 
A\ — + ip 

r 



(101) 



Here, the short distance phase (p is energy independent. 
This can be done by matching the numerical solutions to 
the short distance expanded ones, a cumbersome proce- 
dure in practice [44| . It is far more convenient to use an 
equivalent short distance cut-off method with a bound- 
ary condition. Thus, at the cut-off boundary, r = r c we 
can impose a suitable regularity condition depending on 
the sign of g 2 NN — f 2 ^N- ^ set °f possible auxiliary 
boundary conditions was discussed in Ref. [44| . showing 
that the rate of convergence was depending on the partic- 
ular choice. Actually, there are infinitely many auxiliary 
boundary conditions which converge towards the same 
renormalized value, as we discuss below. 



B. The deuteron 

In this case we have a negative energy state 

,,2 



E = - 1 -. 
M 



(102) 



and we look for normalized solutions of the coupled equa- 
tions (|94[) and normalized to unity, 



dr [u{r) 2 + w(r) 2 ] = 1 , 



(103) 



which asymptotically behave as 



u 7 (r) 



Asve- ir 1 + 



j r (7r) ; 



(104) 



where As is the asymptotic wave function normalization 
and 77 is the asymptotic D/S ratio. To solve this problem 
we introduce, as suggested in [44], the auxiliary problems 




~7r 



3 

7r 



(7& 



(105) 
(106) 



which solutions depend on the deuteron binding energy 
through 7 and the OBE potential. Further, we can use 
the superposition principle of boundary conditions to 
write 



u(r) = u s (r) + r)u D (r) , 
w(r) — ws(r) + rjwoir) . 



(107) 



The short distance regularity conditions (see below) must 
be imposed an a cut-off radius r c in order to determine 
the value of f](r c ). Then, for a given solution we compute 
several properties as a function of the cut-off radius, r c . 
From the normalization condition, Eq. (| 1 03|) , in r c < r < 
00 we get As{r c ). In this paper we also compute the 
matter radius, 

(r 2 \ A f°° 
rl = Y=A] r r2{u{r)2 + w{r)2)dr ' (108) 

the quadrupole moment (without meson exchange cur- 
rents) 

1 f°° 

Qd= 20 r 2 w(r)(2V2u(r) - w(r))dr , (109) 



the D-state probability 

POO 

Pd = I w{r) 2 dr . 



(110) 



which in the impulse aproximation and without meson 
exchange currents can be related to the deuteron mag- 
netic moment. Finally, we also compute the inverse mo- 
ment 

/>oo 

(r~ 1 )= r^iuir) 2 +w(r) 2 )dr, (111) 



which appears, e.g., in the multiple expansion of the ir- 
deuteron scattering length. 

As mentioned, there are infinitely many possible aux- 
iliary conditions. This is an important point which we 
wish to illustrate. For instance, we could take 



sin au(r c ) + cos aio(r c ) = , 



(112) 
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FIG. 8: Short distance cut-off dependence of the asymptotic 
D/S-ratio for the case with tt + a + p + u. We show the 
dependence of the asymptotic D/S normalization rj for several 
choices of the arbitrary and auxiliary short distance condition 
sin au(r c ) + cos aw(r c ) = for several values of a. 



where we may choose the parameter a arbitrarily 16 . 
This is illustrated in Fig. ([5]). Note that despite possible 
wild behaviour all choices converge to the same value, 
although at a quite different rate. This is indeed an- 
other reason for removing the cut-off although it may be 
appealing and less demanding to choose one particular 
scheme where stability is found at the largest possible 
distances. 

Here we will take the smoothest auxiliary condition 
(labeled as BC6 in Ref. Q) 



u'(r c ) - V2w'(r c ) = , g 2 nNN - f 2 pNN > , 
V2u'(r c ) + w\r c ) = , & NN - f pNN < 



(113) 



Clearly, for the values that we will be using the conver- 
gence is determined by the size of the short distance scale 
characterizing the most singular component of the poten- 
tial. As we see from Eq. (|97l) it depends strongly on the 



combination g 2 NN — fp^N- This is an important point 
since the short distance cut-offs, r c , for which conver- 
gence is achieved may change by orders of magnitude . 
An additional numerical problem arises due to undesired 
amplification of the short distance growing exponential, 
setting some limitations to the numerics due to roundoff 



16 This arbitrariness is not exclusive to this boundary condition, it 
is also present when the standard from factor regularization is 
introduced. The exponential, Eq. ( 1311 . and monopole Eq. (130 ft 
form factors are just two possible choices which do not cover the 
most general form which might allow a theoretical estimate on 
the systematic error. 

17 An extreme example is given by the exceptional case f p MN 

'}.. v v since the 1/r 3 singularity turns into a slowly and loga- 
rithmically converging Coulomb singularity. This case is treated 
specifically in Appendix [B] 



errors. In all our calculations we have payed particular 
attention to these delicate issues. 

The cut-off dependence of these observables is shown 
in Fig. [SI for the case of ir only (Ref. [4J]), tt + a and 
7r + a + p + oj and as we see good convergence can be 
achieved as r c — > 0. As already mentioned, the rate of 
convergence depends on the scale of the singularity. 

The resulting coordinate space deuteron wave func- 
tions, u and w, are depicted in Fig. [TU] for the case of 
tt only (Ref. [4J|), n + a and n + a + p + lo and com- 
pared to the wave functions of the high quality Nijmegen 
potential [6j|. As we see, after inclusion of the scalar 
and vector mesons, the agreement is quite remarkable 
in the region above 1.4 — 1.8fm, their declared range of 
validity. Similarly to the singlet case, we observe oscilla- 
tions in the region below 1 fm. The first node is allowed 
since we are dealing with a bound state, the second node 
occurs already below 0.5fm indicating, similarly to the 
channel, the appearance of infinitely many spurious 
bound states, as we see from the short distance oscillatory 
behaviour of the wave function, Eq. (|101|) . To compute 
such states we proceed similarly to the singlet channel. 
We solve Eq. ([9"4"| with negative energy Eb = —j B /M, 
the asymptotic behaviour in Eq. (|104p and impose the 
regularity conditions, Eq (|113|) . as well as orthogonality 
to the deuteron state, namely 



u^u'g — u ' ub + w^w'g 



= 0, (114) 



from which "/b can be determined. For instance, for the 
scalar parameters in Eq. ([84)) and f P NN = 15.5 we iden- 
tify the first spurious bound state {ub\, Wbi) having one 
node less than the deuteron wave functions (ud, Wd) tak- 
ing place at ^bi = 3.438fm~ . The corresponding en- 
ergy is Ebi — —^i"b\/M = — 490MeV, S-wave normal- 
ization Abi = IS.h&hrT 1 / 2 matter radius tbi = 0.49fm 
and asymptotic D/S ratio r\B\ = 0.1656. This state is 
clearly beyond the range of applicability of the present 
framework. Subtracting this pole to the 3 Si amplitude 
would result, according to Eq. lj92|L in Aro = O.Olfm. 
The next spurious state has Eb2 < — 18GeV!!. Note 
that if the scale where the second unphysical node takes 
place was to be interpreted as a ("hard core") break- 
down distance scale of our approach for the deuteron, it 
is certainly beyond the accessible region at the maximal 
energy in elastic NN scattering. This issue is relevant 
for the calculation of phase shifts where such oscillations 
also occur. The variation of the observables from this 
breakdown scale to the origin, could be interpreted as 
a source of systematic error coming from the fact that 
there is only one bound state and not infinitely many. 
As we see from Fig. [9] the effect is indeed small. 

Numerical results for renormalized quantities can be 
looked up in Table [Til As we see, the inclusion of a 
provides some overall improvement while p and oj yield a 
fairly accurate description of the deuteron for the choice 
fpNN = 15-5 and g^NN = 9 (this latter value complies 
to the SU(3) relation g^NN = 3gpNN when g P NN ~ 2.9). 
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FIG. 9: Short distance cut-off dependence of deuteron properties for the cases with tt, n + a and it + a + p + ui. We show the 
dependence of the asymptotic D/S normalization 77 (upper left panel), the S-wave normalization As (in fm _1//2 , upper middle 
panel), the matter radius r m (in fm, upper right panel), the quadrupole moment Qd (in fm 2 , lower left panel), the D-state 
probability (lower middle panel) and the inverse radius (r ) (in fm lower right panel). Experimental or recommended values 
can be traced from Ref. 16711. 
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FIG. 10: Deuteron wave functions, u (left) and w (right), as a function of the distance (in fm) in the OBE. We show lir, ln + la 
and 7r + <7 + p + iu compared to the Nijmegen II wave functions [63|. The asymptotic normalization u — > e~ lr has been adopted 
and the asymptotic D/S ratio is obtained to be rj w — 0.2633 and rj naujp = 0.2597 (see table HI)) . 



We show in Fig. (fTTj) the dependence of several proper- 
ties when both the vector mesons p and uj are simultanc- 
osuly considered. In Fig. QT]we plot the dependendcncc 
of (renormalized) deuteron properties as a function of 
fpN n for several values of the effective coupling constant 

g^NN = JdlNN - fpNN m2 P / 2M N featuring the strong 



correlation in the 1 5'o channel pointed out in Section fVl 
The scalar coupling g a NN and scalar mass m a are al- 
ways readjusted to fit the 1 S'o phase shift. As we see, 
for the asymptotic D/S ratio, 77, there is a wide range 
of possible values within the experimental uncertainties 
but we obtain the bounds J p nn < 15 and g^NN < 15- 
It is amazing that the value of the tensor-p coupling is 
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FIG. 11: Dependence of the deuteron observables as a function of / p jvjv for several values of the the effective coupling constant 

OojNN — \JotNN ~ fpNN m p/'2M%. ga-NN and m a are always readjusted to fit the 1 So phase shift. We show the dependence 

of the asymptotic D/S normalization 77 (upper left panel), the S-wave normalization As (in fm -1 / 2 , upper middle panel), the 
matter radius r m (in fm, upper right panel), the quadrupole moment Qd (in fm 2 , lower left panel), the D-state probability 
(lower middle panel) and the inverse radius (r _1 ) (in fm -1 lower right panel). The leading N c contributions to the OBE 
(ct + 7r + p + oj) potential are considered. Experimental or recommended values can be traced from Ref. [oTj . 



TABLE II: Deuteron properties and low energy parameters in the 3 Si — 3 D 1 channel for OBE potentials including ir, tv + <r, 
7r + a + p + ui. We use the non-relativistic relation 7 = ^/2p nv B with B = 2.224575(9) and take m = 138.03MeV, and 
g-wNN = 13.1083 [gj. From a fit to the 1 5'o channel we have m CT = 501MeV and g a NN = 9.1. The simplifying relation 
gujNN = f P NNrn p /\f2MN is used throughout, irapuj corresponds to take / p jvjv = 15.5 and g^NN = 9.857 while napui* 
corresponds to take / p nn = 17.0 and g^NN = 10.147 
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"(Non relativistic). See e.g. Ref. (6711 and references therein. 



so well determined to be f P NN ~ 16 — 17 and corre- 
sponds to the strong k p situation described by Machleidt 
and Brown [68] . Note that results depend in a moderate 
fashion on f P NN for not too large values, as one would 
expect from the short range of the p— meson. 



C. Zero energy 

At zero energy, the asymptotic solutions to the coupled 
equations (|94|) are given by 
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FIG. 12: Dependendence of the low energy parameters as a 
function of / p jvjv for several values of the the effective cou- 
pling constant g^NN 



\J91nn - fpNN m2 P / 2M N- 9°nn and 
ma are always readjusted to fit the So phase shift. We show 
the dependence of the 3 Si scattering length ao (in fm) and 
effective range ro (in fm) as well as Q02 (in fm 3 ) and 0:2 (in 
fm ) compared to the experimental values or the Nijm2 and 
Reid93 potentials (horizontal straight lines) The leading JV C 
contributions to the OBE (a + n + p + u) potential are con- 
sidered. 
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where ao, «2 and ao2 are low energy parameters obtained 
from the phase shifts (see Sect. I VI Dp . Using these zero 
energy solutions one can determine the effective range. 
The 3 5i effective range parameter is given by 



r = 2 
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u , a (r) 2 - w . a (rY 



dr 



(116) 



Using the superposition principle of boundary conditions 
we may write the solutions as 



uo, a (r) 
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where the functions ui,2, 3, 4 and wi.2,3.4 are independent 
on ao, «02 and a 2 and fulfill suitable boundary condi- 
tions. The orthogonality constraints for the a and (3 
states read in this case 



= 

= 



(118) 



A further condition which should be satisfied is the a — (3 
orthogonality 



UQ.aU /3 - M q Mo,/3 + Wo, a W /3 - W a UQ,f3 



= 

(119) 



as well as the short distance regularity conditions, 
Eq. (|113[) . In all we have an over-determined system with 
5 equations and three unknowns «02, &2 and ao- Solving 
the equations in triplets we have checked the numerical 
compatibility at the 0.01% level for the shortest cut-offs, 
r c ~ 0.02fm typically used. Using the superposition prin- 
ciple decomposition of the bound state, Eq. (|107p . and 
for the zero energy states, Eq. (|117p . one can make the 
orthogonality relations explicit in ao, ao2, ot 2 . The val- 
ues of CV02 and a 2 are not so well known although they 
have been determined from potential models in Ref . [69j . 

In Fig. [T2] we show the dependence of the low en- 
ergy parameters of the leading N c contributions to the 
OBE (cr + 7r + p + lu) potential as a function of / p nn 
for several values of the the effective coupling constant 
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always readjusted to fit the 1 S'o phase shift. Similarly 
to the deuteron case we observe stronger dependence on 
fpNN and a relative insensitivity on the effective cou- 
pling g^NN- We remind that along any of these curves 
the 1 So phase shift is well reproduced with an acceptable 
X 2 /DOF < 1. As we see, the values J p nn = 17.0 and 
9uNN = reproduce quite well the low energy parame- 
ters, corresponding to the reasonable g U NN — 10.4. 

Numerical results for the low energy parameters are 
shown in Table [TTJ Again, the inclusion of a provides 
some overall improvement while p and lu yield a better 
description of the deuteron for the choice J p nn = 15.5 
and gu>NN — 9.0. There is nonetheless a small mismatch 
to the experimental or recommended potential values 
when the zero energy wave functions are obtained from 
the orthogonality relations to the deuteron, Eq. (|118p . 
u i{ r ) 7 As one can see further improvement is obtained when 
fpNN = 17.0 and g^NN — 10.3. In this case we get a 
w^{r) , SU(3) violation; g^NN = ^-^9pNN, which actually agrees 
with the expectations from radiative decays lu — > e + e~ 
(117)and p — > e + e~ (see e.g. [70l|). 
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FIG. 13: Convergence of the np Spin Triplet Eigen phase shifts for the total angular momentum j — 1 as a function of the 
short distance cut-off radius r c (in fm) for several fixed values of the CM momentum p — 100, 200, 300 and 400MeV. 




FIG. 14: np Spin Triplet Eigen phase shifts for the total angular momentum j = 1 as a function of as a function of the CM 
momentum. We show tt, tt + a and tt + a + p + uj compared to an average of the Nijmegen partial wave analysis and high 
quality potential models [H], HU]. We take (/ p jvjv, Q^nn) = (15.5,9.857), ( fp n n , Quj nn} = (17.0, 10.147), . 



D. Phase shifts 

Finally, in the case of positive energy we consider 
Eq. JMD with 



E 



El 

M ' 



(120) 



with p the corresponding CM momentum. We solve 
Eq. ([94]) for the a and (3 positive energy scattering states 
and choose the asymptotic normalization 



Uk, a {r) 
w k , a {r) 
u k,f3(r) 



sin S 



-Uo(kr) 
l v 



cosJi — yo(kr) sin 5] 
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anii \ 

- ■ , (jo (At) cos 8 2 - yo (kr) sin 6 2 
sin 01 V 

-^r(h(kr) cos<5 2 - faikr) sin^ 2 
iindi V 



(121) 



where ji (x) = xji (x) and yi (x) — xyi (x) are the reduced 
spherical Bessel functions and £1 and 82 are the eigen- 
phases in the 3 S\ and 3 D\ channels, and e is the mixing 
angle E\. To carry out the renormalization program, we 
use the superposition principle of boundary conditions 
which makes the discussion more transparent. Let us 
define the four auxiliary problems 
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(122) 



which depend solely on the potential and can be obtained by integrating in. Thus, the general solution satisfying the 
a and (3 asymptotic conditions can be written as 

4 4 4 4 

Uk,a(r) = c i<a u k ,i(r) , w k<a (r) =^2ci <a w kt i(r) , u k ,p{r) = ^ <k,pUk,i{r) , w k ,p(r) = ^ Cj t pw kti (r) ■ (123) 



Fixing the constants to the asymptotic conditions 
Eq. {TUP we get 
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(124) 

In the low energy limit e — ► — ao2& 3 5 a - ► —ot^k and 
5/3 — > — (a2 — a^/ao)^ 5 and the zero energy solutions 
discussed in Sec. IVICI are reproduced. The use of the 
superposition principle for boundary conditions as well 
as the orthogonality constraints to the deuteron wave 
analogous to Eq. (|72f yields 







U l u k,a ~ « 7 «*,a + W 7 W k a - W 7 U k , a 

u i u 'k,p - u' 7 u kt p + w^w' h p - w'^u^p 



= 



(125) 



which together with the short distance regularity con- 
ditions, Eq. (|113[) allow us to deduce the corresponding 
3 Si — 3 D\ phase-shifts. A further condition is the a — (3 
orthogonality 



U k , a U k p - v! Ka U k> p + Uli, a Wy - W k ,a U k,P 



= 
(126) 



In all we have an over-determined system with 5 equa- 
tions and three unknowns. We have checked that almost 



any choice yields equivalent results with an accuracy of 
0.001° for the highest CM momenta and the shortest cut- 
off, r c ~ 0.02fm. 



As we have mentioned already, the numerical solution 
of the problem requires taking care of spurious amplifi- 
cation of the undesired growing exponential at any step 
of the calculation. The situation is aggravated by the 
fact that for the phase shifts the maximum momentum 
p = 400MeV explores the region around 0.1 — 0.5fm, so 
it is important to make sure that we do not see cut-off 
effects in this region. To provide a handle on the nu- 
merical uncertainties we show in Fig. (TIB")) the results 
for the phase shifts Si, 82 and e as a function of the 
cut-off radius, r c and for several fixed CM pn momenta, 
p = 100, 200, 300, 400MeV. As we see, there appear clear 
plateaus between 0.1 — 0.2fm which somewhat steadily 
shrink when the momentum is increased. Note that 
these values of the short distance cut-off translates into 
a CM momentum space cut-off range A = 7r/(2r c ) = 
1.5-3GeV. 



The results for the 3 Si — 3 Di phase shifts as a function 
of the CM momentum are depicted in Fig.[JJ]for tt, it + <t 
and 7r + a + p + oj and compared to the Nijmcgcn anal- 
ysis 0, HI]. We use g aNN = 9.1, ra a = 501MeV and 
when vector mesons are included we take J p nn = 15.5 
and (JujN n = 9 or f pNN = 17.0 and g^NN = 10.147 
corresponding to Sets irapui and irapuj* in Table [TT] re- 
spectively. On a first sight we see an obvious improve- 
ment in both the 3 Si and 3 Di phases and not so much 
in the mixing angle E\ as compared to the simple OPE 
case. One should note, however, that besides describing 
by construction the single phase shift 1 5'o (see Fig. (0| 
) we also improve on the deuteron (see Table HI)) . Obvi- 
ously, it would be possible to provide a better description 
of triplet phase shifts, however, at the expense of wors- 
ening the deuteron properties and the singlet channel. 
Clearly, there is room for improvement, and our results 
call for consideration of sub-leading large N c corrections 
in the OBE potential. This would incorporate, the rel- 
ative to leading l/N% relativistic corrections, spin-orbit 
effects, finite meson widths, non-localities, etc. 



23 




FIG. 15: Convergence of the np Spin Singlet 1 S'o phase shift (in degrees) as a function of the short distance cut-off radius r c 
(in fm) for several fixed values of the LAB energy when the scattering length ag = — 23.74fm is also fixed. We compare the 
purely renormalized case with no form factors (solid line) with the renormalization also including the exponential (dotted line) 
and monopole (dashed line) form factors and for the cut-off values A^jvat = 1300MeV (left panel) and A^-atjv = 2000MeV (right 
panel), all others fixed to A&nn = A p jvjv = A^nn ~ 2GeV. We also add the error bands related to the Nijmegen PWA and 
high quality potentials |l4, 63]. All other meson parameters in the OBE potential are kept the same. 



TABLE III: Fits to the Renormalized 1 Sq phase shift of the Nijmegen group using the OBE potential without or with 
strong exponential and monopole form factor. We fix ao — — 23.74fm and take m = 138.03MeV, and g^NN = 13.1083 [64[ and 
m p — = 770MeV and fit m CT , g a NN and fix g^ NN = 0. We use A^jviv = 1300MeV and A ct jvjv = A p nn = A^nn = 2000MeV. 
Eb represents the energy of the (spurious) bound state when it does exist. 
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VII. INFLUENCE OF STRONG FORM 
FACTORS IN THE RENORMALIZATION 
PROCESS 

Given the reasonable phenomenological success of the 
renormalization approach one may naturally wonder 
what would be the effect of the form factors in our calcu- 
lation. In this section we discuss the influence of strong 
form factors in the calculated properties on top of the 
renormalization process. Our main quest is to find out 
whether they lead to observable physical effects after 
renormalization. An equivalent way of posing the ques- 
tion is to determine whether finite nucleon size effects can 
be disentangled explicitly in NN scattering in the elastic 
region. 

To analyze this important issue in detail, in Fig. [15] 
we show the phase shift in the channel for fixed 
LAB energy values as a function of the short distance 



cut-off radius r c when the scattering length is fixed to 
its experimental value, ao = — 23.74fm as we explained 
in Section [V] We use the same parameters as for the 
renormalized solution without vertex function, for sev- 
eral fixed values of the LAB energy and for the cut- 
off values A wA rjv = 1300MeV and A nNN = 2000MeV, 
all others fixed to A ct ;vjv = A p nn = A^nn = 2GeV. 
As one clearly sees strong form factors are invisible for 
r c > 0.3fm. For lower values of the short distance cut- 
off r c both monopole and exponential form factors agree 
with each other but deviate strongly from the Nijmegen 
database. Note that the lines should be supplemented 
with estimates of theoretical errors, not shown to avoid 
cluttering of the plot. When those errors are included 
the Nijmegen data are basically compatible with the the- 
oretical curves in the flat preasymptotic region around 
0.3 — 0.5fm (see also the discussion around Fig. [S]). 

Of course, one may attribute the discrepancy to the 
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FIG. 16: Short distance cut-off r c dependence of deuteron properties for the irapui case (see table [TTJ) . We compare the 
purely renormalized calculation with the cases for both exponential, Eq. (|31|l . and monopole Eq. (|3U|I form factors taking 
A„-]vjv = 1300MeV, all other cut-offs being kept to Ao-tvjv = A p jvjv = A^jviv = 2000MeV. We show the dependence of the 
asymptotic D/S normalization 77 (upper left panel), the S-wave normalization As (in fm -1 / 2 , upper middle panel), the matter 
radius r m (in fm, upper right panel), the quadrupole moment Qd (in fm 2 , lower left panel), the D-state probability (lower 
middle panel) and the inverse radius {r~ } (in fm -1 lower right panel). Experimental or recommended values can be traced 
from Ref. |E3]. 



choice of parameters, which have been chosen to fit the 
renormalized solution without form factors. A some- 
what complementary way of seeing this is by refitting 
the parameters using both exponential and monopole 
vertex functions but fixing by construction the scatter- 
ing length uq. The results for A^yyyy = 1.3GeV and 
Ao-jvjv = A p jva7 = Apjvjv = 2GeV are displayed in Ta- 
ble [TTTJ As we see, the parameters change almost within 
the uncertainties, showing the marginal effect of the ver- 
tex functions after renormalization. Due to the presence 
of non-linear correlations, difficult to handle by standard 
means, we have fixed g^NN to its minimum value (com- 
patible with zero) and estimated its error by varying it 
independently from its mean value to values still giv- 
ing an acceptable fit, yielding g* NN — 0(3). We also 
show the effect of the short distance cut-off r c which, as 
we see, is rather small. Overall, these results provide a 
further confirmation of our naive expectations; nucleon 
finite size effects and vector mesons do not provide the 
bulk in NN scattering in central waves, and actually can- 
not be clearly resolved. Of course, this should be checked 
in higher partial waves, but those are expected in fact to 
be less sensitive to short distances. 



Finally, in Fig. [TH] the influence of the vertex functions 
is analyzed for some of the computed deuteron properties. 
As we see there is a fair coincidence of the purely renor- 
malized solution with no form factors with the equally 
renormalized solution including the form factors in the 
potential in the region around r c ~ 0.3 — 0.6fm. The 
deviation below 0.3 fm signals the onset of the irregular 
D-wave solution, which behaves as w(r) ~ r~ 2 at small 
distances and hence yields eventually a divergent result. 
Note that in order to have a smooth behaviour at short 
distances when renormalization is over-imposed to the 
potential with form factors we should choose the regular 
D-wave solution w(r) ~ r 3 but then the potential param- 
eters, either couplings or form factor cut-off parameters 
should also be fine-tuned. 

While it is fairly clear that vertex functions do ex- 
ist and are of fundamental importance, it is also true 
that they start playing a role as soon as the probing 
wavelength resolves the finite nucleon size. Our calcu- 
lations suggest on a quantitative level that provided the 
NN scattering data are properly described with form fac- 
tors, they will be effectively irrelevant under the renor- 
malization process, and for CM momenta below 400MeV, 
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vertex functions are expected to play a marginal role. 



VIII. CONCLUSIONS AND OUTLOOK 

In the present paper we have analyzed the OBE po- 
tential from a renormalization point of view. As we 
have shown, the meson-nucleon Lagrangean does not pre- 
dict the S-matrix beyond perturbation theory. The non- 
perturbative nature of low partial waves and the deuteron 
in the NN problem suggests resuming OBE diagrams by 
extracting the corresponding potential. The OBE poten- 
tial, however, presents short distance divergences which 
make the solution of the corresponding Schrodinger equa- 
tion ambiguous. The traditional remedy for this problem 
has been the inclusion of phenomenological form factors 
which parameterize the vertex functions in the meson 
exchange picture. We have shown that the meson ex- 
change potential with form factors generates spurious 
deeply bound states for natural values of the coupling 
constants. The price to remove those is to fine tune the 
potential at all distances, and in particular at short dis- 
tances. Thus, while it is claimed that vertex functions 
implement the finite nucleon size, it is very difficult to 
disentangle this from meson dressing and many other ef- 
fects where the meson theory does not hold. 

The renormalization approach suggests that extracting 
this detailed short distance information may in fact be 
unnecessary for the purposes of Nuclear Physics and the 
verification of the meson exchange picture. Contrarily to 
what one might naively think, renormalization is a prac- 
tical and feasible way of minimizing short distance am- 
biguities, by imposing conditions which are fixed by low 
energy data independently on the potential. We have ar- 
gued that within this approach we face from the start our 
inability to pin down the short distance physics below the 
smallest de Broglie wave length probed in NN scattering. 
Indeed, the central scattering waves and the deuteron can 
be described reasonably well and with natural values of 
the meson-nucleon couplings. Within the standard ap- 
proach this could only be achieved by fine tuning meson 
parameters or postulating the meson exchange picture 
to even shorter ranges than 0.5fm. In our case the inclu- 
sion of shorter range mesons induces moderate changes, 
due to the expected short distance insensitivity embodied 
by renormalization, despite the short distance singularity 
and without introducing strong meson-nucleon-nucleon 
vertex functions. If phenomenological vertex functions 
are added on top of the renormalized calculations mi- 
nor effects are observed confirming the naive expectation 
that finite nucleon size ~ 0.5fm need not be explicitly in- 
troduced within the OBE calculations for CM momenta 
corresponding to the minimal wavelength l/p ~ 0.5fm. 

The renormalization process introduces spurious 
deeply bound states regardless on whether or not the 
potential is regular or singular. This can be appreciated 
in the excessive number of nodes of the wave function 
close to the origin, in the region below 0.5fm. We have 



checked that the corresponding CM energies are in ab- 
solute values much higher than the maximum scattering 
CM energies, and hence the role played by these spurious 
states is completely irrelevant. We remind that within 
the standard approach with form factors those spurious 
bound states also take place when natural values of the 
coupling constants are taken. 

One of the problems with potential model calculations 
is the ambiguity in form of the potential, since it is de- 
termined from the on-shell S-matrix in the Born approxi- 
mation and an off-shell extrapolation becomes absolutely 
necessary. In the large N c limit the spin-isospin and kine- 
matic structure of the NN potential simplifies tremen- 
dously yielding a non-relativistic and uniquely defined 
local and energy independent function. Relativistic ef- 
fects, spin-orbit, non-localities as well as meson widths or 
other mesons enter as sub-leading corrections to the po- 
tential with a relative order 1/AT?. However, it consists of 
an infinite tower of multi-meson exchanged states, which 
range is given by the Compton wavelength of the total 
multi-meson mass. One of the advantages of the large 
N c expansion is that it is not particularly restricted for 
low energies. This is exemplified by several recent calcu- 
lations of NN potentials using the holographic principle 
based on the AdS/CFT correspondence [Hi, [72, III 18 ■ 
A truncation of the infinite number and range of ex- 
changed mesons is based on the assumption that the 
hardly accessible high mass states are irrelevant for NN 
energies below the inelastic pion production threshold. 
This need not be the case, unless a proper renormaliza- 
tion scheme makes this short distance insensitivity mani- 
fest. Actually, within such a scheme the counterterms in- 
clude all unknown short distance effects, but enter as free 
parameters which do not follow from the potential and 
which must be fixed directly from NN scattering data or 
deuteron properties. In the present work we have imple- 
mented a boundary condition regularization and carried 
out the necessary renormalization. This allows, within 
the OBE potential to keep only n, cr, p and u> mesons 
and neglect effectively higher mass effects for the lowest 
central s- waves as well as the deuteron wave function. In 
many regards we see improvements which come with very 
natural choices of the couplings, and are compatible with 
determinations from other sources. From this viewpoint, 
the leading N c contribution to the OBE potential where 
7r, cr, p and u> mesons appear on equal footing, seems su- 
perior than the leading chiral contribution which consists 
just on 7r. 

The value of the a mass was fixed by a fit to the 
1 S phase shift yielding m a = 501(25)MeV. The val- 
ues obtained from the coupling constants reproducing 
the 1 5'o and 3 Si — 3 D\ channels are very reasonable tak- 



In this calculations only tt,p,lj and A\ mesons and their radial 
excitations contribute. Note, however, that the only contribution 
to the central force Vc stems from the tower u,u',u", . . . which 
is generally repulsive. 
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ing into account the approximate nature of our calcula- 
tion, g aNN = 9(1), QuN n = 9.5(5) and f pNN = 16.3(7); 
the range is compatible with the putative 10% accu- 
racy of the 1/N% corrections. For the accepted value 
9 p nn = 2.9(1) this yields guiNN / 9pNN = 3.27(17) a value 
in between the SU(3) prediction g^NN / 9pNn\su(3) — 3 
and the one from the e + e~ — * p and e + e~ — * uj decay ra- 
tios, QuiNN / 9pN n\c+ e~ = 3.5. We also get Jpnn/ 9pNN = 
Hp = 5.6(3); a value in agreement from tensor coupling 
studies. It is noteworthy that the repulsion triggered 
by the to meson is not as strong and important as re- 
quired in the conventional OBE approach where usually 
a strong violation of the SU(S) relation is observed as 
well. The reason is that, unlike the traditional approach, 
the renormalization viewpoint stresses the irrelevance of 
small distances. This is done by the introduction of coun- 
terterms which are fixed by threshold scattering param- 
eters at any given short distance cut-off scale r c . For the 
minimal de Broglie wavelength probed in NN scattering 
below pion production threshold, 1/p ~ 0.5fm, a stable 
result is obtained generally when r c = 0.1 — 0.2fm. Any 
mismatch to the observables can then be attributed to 
missing physical effects. While the present calculations 
are encouraging there is of course room for improvement. 

One serious source of complications and limitations for 
renormalization in general lies in its difficult marriage 
with the variational principle [74j . The existence of two- 
body spurious deeply bound states drives naturally the 
energy of the system to its lowest energy state, if allowed 
to. On the other hand, one should recognize that the ex- 
istence of a minimum is tightly linked to a subtle balance 
between kinetic and potential energy, which undoubtely 
exists but may well take place beyond the applicability 
range of the meson exchange picture requiring an artifi- 
cial fine tuning. This clearly influences the three, four, 
etc. body problems if they would be treated in the stan- 
dard and variational fashion but not necessarily so if the 
few body problem is consistently renormalized. Our re- 
sults show that one has to choose between fine tuning 
and renormalization. The standard approach has tradi- 
tionally been sensitive to short distance details and has 
required fine tuning meson coupling constants, in partic- 
ular those corresponding to vector mesons, to unnatural 
values. In contrast, the renormalization approach is free 
of fine tuning, and allows to fix meson constant from 
other sources to their natural values. 

While we have been using the leading large N c con- 
tributions to the full OBE as a simplifying book-keeping 
reduction, we do not expect that such an approximation 
becomes crucial regarding the main conclusions on form 
factors. However, the most speculative prospective of 
the present calculation lies in the possibility of promot- 
ing it to a model independent large N c result. One should 
bear in mind, however, that we have only kept leading 
N c OBE contributions. There is, of course, the delicate 
question on what 2tt, 3tt and A contributions should be 
considered, firstly to avoid double counting with the col- 
lective a, p and w states, and second to comply with 



the large N c requirements. To our knowledge, the expec- 
tations of Ref. [53] of a large A c consistent multimeson 
exchange picture have not been explicitly realized for the 
chiral potentials without [75[ and with [76| A-isobar con- 
tributions as they do not scale properly with N c ; one 
has 9a ~ N c , /tt ~ y/N c and there are terms scaling 
as V r 2 < ^ hPT ~ g\j ~ N% and not as ~ A c as found in 
Refs. [i^, [52|] ■ Our results suggest a scenario where the 
multimeson contributions invoked in Ref. [53| would in- 
deed be small, but this should be checked explicitly. One 
further complication comes from the fact that in the large 
N c limit the nucleon-delta splitting becomes small, and in 
fact lighter than the pion mass. According to the Regge 
theory formula M\ — Mjy = m 2 p — m 2 [77| and assuming 
the scaling M/v = N c m p /2, the crossover between both 
mass parameters happens at about N c ~ 6. Actually, in 
the strict limit one should consider not only NN but at 
least also NA and AA channels as well, as they become 
degenerate. The calculation of [52|, UM 0ln y includes the 
restriction of the baryon-baryon interaction to the NN 
sector. In a more elaborate treatment one should include 
the A as intermediate dynamical states which in the elas- 
tic NN region contribute as sub-threshold effects [79| 
which decouple for large AA splitting but which become 
degenerated when the A A splitting is driven to zero. 
In addition, it would also be interesting, still within the 
OBE framework, to see what is the effect of the relative 
1/A^ corrections, which include in particular relativis- 
tic, non-local, finite-size, spin-orbit, finite meson width 
corrections as well as other mesons. 

Finally, let us also note that besides the many im- 
provements mentioned above to the present calculation, 
the possibility of making a good phenomenology while 
avoiding strong form factors in the NN potential has 
further and important benefits. In particular it makes 
the discussion of gauge invariance much simpler, as we 
are effectively dealing with local theories with no cut- 
off. Under this circumstance the cumbersome gauging 
procedures involving path-dependent link operators and 
which becomes necessary in order to minimally imple- 
ment gauge invariance would not be needed. In a recent 
communication [571 ] we have evaluated electromagnetic 
deuteron form factors in the impulse approximation and 
using the renormalization scheme presented in this pa- 
per, with a reasonable momentum transfer dependent 
behaviour up to about q ~ 800MeV and definitely im- 
proving over OPE. Actually, these form factors as well 
as some of the presently computed deuteron properties 
are expected to have significant corrections from MEC. 
Let us remind that MEC are a genuine consequence of 
the Meson Exchange picture in the NN interaction, but 
in fairness also require constructing exact NN wave func- 
tions from the corresponding Hamiltonian, as we have 
done here. The present paper shows that renormalization 
for the OBE potential is not only feasible as a previous 
and theoretically appealing step to evaluate matrix ele- 
ments of electroweak currents but also and perhaps sur- 
prisingly yields a sound phenomenologically. It also helps 
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in reducing the impact of the hardly accesible short dis- 
tance region of the nucleon-nucleon interaction, thereby 
reducing standard and much debated ambiguities. It re- 
mains to be seen if this holds true also for low energy 
electroweak reactions where the meson exchange picture 
is traditionally expected to work. 
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APPENDIX A: OVERVIEW OF COUPLING 
CONSTANTS 

A crucial point in the present framework corresponds 
to the choice of coupling constants, g^NN, QaNN, foNN 
and g^NN (for an older review see e.g. Ref. [7(j]) and 
masses, m 7!l m CT , m p and m u , entering the calculation. 
We review here reasonable ranges on the basis of several 
sources, but bearing in mind that we are keeping only 
the leading N c contributions to the OBE potential. 

• 9-kNN- According to the Goldberger-Treiman re- 
lation (subjected to pion mass corrections and/or 
higher meson states) the pion nucleon coupling 
constant should be g-^NN — QaMn/S-k = 12-8 
for the axial coupling constant gA = 1.26. A 
phase shift analysis of NN scattering [64( yields 
g-KNN = 13.1083. Nevertheless, the latest deter- 
minations from the Goldberger-Miyazawa-Ochmc 
(GMO) sum rule [8(| yields the value g^NN = 
13.3158 this variation at the 5% level dominates 
the uncertainties in the Itt exchange calculations. 

• gaNN- For the scalar coupling constant, the 
Goldberger-Treiman relation for scalar mesons 
yields, g a NN — Mn/ f-n — 10.1. However, if we con- 
sider contributions from excited scalar mesons we 
may expect a a somewhat different number. Actu- 
ally, QCD sum rules yield g aNN = 14.4±3.7 for 
the Ioffe current nucleon interpolator and a smaller 
value g^NN = 7 ± 3 for more general interpola- 
tors [821 . A rece nt q uark model calculation yields 
g aN N = 14.5 ± 2 df. 

• gpNN- The vector g P NN coupling constant is af- 
ter Sakurai's universality g P NN — gpirir/2 while 
the current-algebra KSFR relation provides <? p7r7r = 
m p /(V2U), yielding g pNN = 2.9. The pNN Ver- 
tex in Vector Dominance Models was also deter- 



mined in the old analysis [84[ yielding yields g p = 
2.9(1) a value confirmed in Ref. (85| . 

• fpNN- The tensor f P NN coupling is usually given 
by the ratio to the vector coupling/^ n = n p g p NN- 
In single vector meson dominance models k p = /i p — 
fjb n — 1 with /i p = 2.79 and fx n = —1.91 the magnetic 
moments (in nuclear magneton units e/(2M p )) of 
proton and neutron respectively, yielding n p = 3.7 
and hence f P NN = 10.7(4) for g P NN = 2.9(1). 

• gu>NN- The relation g U NN — 3g P NN (= 8.7(3) for 
g P NN — 2.9(1)) is the SU(3) prediction for the 
ideal w — <f> mixing case corresponding to the OZI 
rule where g^NN = as well. Vector meson e.m. 
decays u> — > e + e _ and p — > e + e~. account for 
SU(3) breaking as g^NN = 3-5g P NN (= 10.2(4) for 
g pNN = 2.9(1)). 

• fu:NN- The tensor f^NN coupling is also given 
by its the ratio to the vector coupling/^ATjv = 
KujguiNN- In single vector meson dominance models 
Kui — + (J-n — 1 yielding k u = —0.12 and hence 
Unn = -0.3(1) for g uN N = 3 - 3.5. 

Nucleon Electromagnetic Form factors with high en- 
ergy QCD constraints alsoprovide information on vec- 
tor meson couplings. Ref. [86| yields g U NN — 20.86(25) 
and fuiNN = —3.41(24) and k p — 6.1(2), and more re- 
cently 87] it was found g U NN = 20(3) and f^NN = 3(7). 
On the other hand, QCD sum rules yield for the pNN 
coupling a spread of values is g P NN = 2.4 ± 0.6 and 
f pNN = 7.7 ± 1.9 HH andopArjv = 3.2 ± 0.9 and 
fpNN + g P NN = 36.8 ± 13.0 [If 

Phase-shift analyzes of NN scattering below 160 MeV 
based on the t\ mixing angle were argued to be an 
indication for a stron g te nsor force [90| . an issue fur- 
ther qualified in Ref. [68|]. The strong tensor coupling 
is K p = fpNN / gpNN — 6.1(6) and the weak is k p = 
lip — 1 — |U n = 3.7 corresponding to vector meson domi- 
nance saturated with a single state. Note that the value 
f P NN — g-nNN = 13.1 for which the tensor force 1/r 3 sin- 
gularity disappears corresponds to k p — 4.5(2) a value in 
between weak and strong. 



APPENDIX B: THE EXCEPTIONAL 
NON-SINGULAR CASE 



As already mentioned in Section [TlB] there is an excep- 
tional situation f P NN = g-nNN where the OBE potential 
is not singular, Eq. (fl6|) . and the use of form factors 
would not be necessary. If we keep g^NN = 13.1 that 
means f P NN — 13.1, a not completely unrealistic value 
lying in between the single vector meson dominance es- 
timate and the usual OBE value (see Appendix [A"|l . so 
it is worth analyzing this case separately. Since the sin- 
gularity affects mainly the coupled spin triplet channel, 
one may wonder what would be the consequences for the 
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TABLE IV: Deuteron properties for the exceptional case J p nn = QtvNN of non-singular large N c OBE potentials. In all cases 
we take r c = O.OOlfm. We compare renormalized vs. regular solutions for similar choices of parameters. We use 7 = ^/2/i np Ba 
with B d = 2.224575(9) and take g^ NN = 13.1083, = 138.03MeV, m p = = 782MeV. The fit to the 1 So phase shift gives 
m CT = 501MeV and g^NN = 9.1. Experimental or recommended values can be traced from Ref. [67t ] . 
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deuteron. We will show that our conclusions are not ruled 
out by this exceptional case 19 . 

Note that within the renormalization approach this 
particular situation has been scanned through in Fig. [IT] 
where nothing particularly noticeable happens. Actually, 
at short distances we have a coupled channel Coulomb 
problem where the short distance behaviour can gener- 
ally be written as a linear admixture or regular and ir- 
regular solutions, 



u(r) ~ air + a 2 
w{r) ~ bir 3 + b 2 r~ 2 



(Bl) 



In order to get a normalizable wave function we must 
impose the regular solution for the D— wave, meaning 
62 = 0. The renormalized solution corresponds then to 
fix the deuteron binding energy as explained in detail 
in Section I VII and integrate in with the result that the 
S— wave may have an admixture of the irregular solution. 
The regular solution takes the value a 2 = 0. The bound 
state properties are now predicted completely from the 
potential. 

In practice we deal with arbitrarily small but finite cut- 
offs, r c — > 0. In this situation it is simplest to use the 
superposition principle of boundary conditions given by 
Eq. (|107J) for a given energy or 7. From the regularity 
condition of the D-wave we get 



r c ^H = 3, 



(regular D — wave) 



w(r c ) 

which yields the asymptotic D / S-vdXio 

-3wg(r c ) + r c w' s {r c ) 



2>w D {r c ) + r c w' D (r c ) 



(B2) 



(B3) 



19 A compelling scenario where the singularity cancels might hap- 
pen for an infinite tower of exchanged mesons fulfilling the sum 
rule &NN + dl'NN + ■■■ = fpNN + fp'NN +•••=■ Even if this 
was the case the implications after renormalization are meager. 



This provides a relation between 7 and 77. The renormal- 
ized condition yields an arbitrary value of u at the origin, 
so the energy may be fixed arbitrarily, and thus 



u'(r c ) 



7^ 1, (irregular S — wave) (B4) 



u{r c ) 

The regular solution corresponds to 
u'(r c ) 



u(r c ) 



1, (regular S — wave) 



(B5) 



which in general will not be satisfied by the physical 
deuteron binding energy. Thus, for the regular solution 
we will have either a wrong value of the energy or the 
potential parameters must be readjusted. A value of 
r c = O.OOlfm proves more than enough. 

Numerical results for a fixed parameter choice with 
9uNN = are presented in table [IV] As we see the regu- 
lar solution generates a bound state with Eb ~ — 16MeV 
which is clearly off the deuteron with equally bad prop- 
erties. In order to achieve the correct deuteron binding 



energy we just increase the coupling to g* 



iVJV 



3.75 in 



the regular solution case. In this case both renormalized 
and regular solution would coincide accidentally. How- 
ever, if we increase to twice this value g^ NN — 2 x 3.75 
we observe tiny changes in the deuteron properties as 
compared to the g^NN = case when the renormalized 
solution is considered whereas the regular solution be- 
comes unbound. These results illustrate further the sharp 
distinction between regular and renormalized solutions 
where one chooses between fine tuning and short dis- 
tance insensitivity respectively. The corresponding wave 
functions to both the renormalized and regular solutions 
with the same meson parameters are depicted in Fig. 1171 
In both cases inner nodes of the wave functions exihibit 
the existence of deeply bound states, as dictated by the 
oscillation theorem. 

Finally, we might try to analyze the consequences of 
taking V^5 1 (r) = Vis (r) in the exceptional case f P NN — 
9nNN = 13.1 and other parameters from the case with 
no form factor, T = 1, of table [T] for the 1 5'o channel. 
Let us remind that two possible scenarious arise in such 
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Non-singular (f pNN = g„ NN ) 
deuteron wave function u(r) 



Nijmll 

Regular (g* roNN = 0) 
Renorm (g* roNN = 0) 



0.5 



1 



1.5 
r [fml 



2.5 





0.8 




0.6 




0.4 




0.2 













-0.2 




-0.4 



0.5 



Nijmll 

Regular (g*^ = 0) 
Renorm (g* fflNN = 0) 




Non-singular (f pNN = g^) 
deuteron wave function w(r) 



1.5 
r[fm] 



2.5 



FIG. 17: Normalized Deuteron wave functions, u (left) and w (right), as a function of the distance (in fm) in the OBE for the 
exceptional non-singular case / p nn = g-rrNN- We show n + a + p + ui both renormalized and the regular solution with the same 
parameters g^ NN = 0. We compare to the Nijmegen II wave functions [63j(see table ITV) ). 



a case, one with no boud state and another one with a 
spurious deeply bound state. For the 3 Si — 3 Di channel, 
this complies to the standard picture that the deuteron 
becomes bound due to the additional binding introduced 
by the small tensor force mixing with the D— wave, ba- 
sically shifting the S— wave potential to an effective one 
Vs Sl {r) ~ Vi So {r) + W T (r) 2 /V SDl (r) . While in the case 
with no spurious bound state for the 1 Sq we do not get 
any deuteron bound state, in the case with the spurious 



bound state the binding energy is Eb ~ — 50McV. This 
is another manifestation of the fine-tuning discussed at 
length in Section fill Bl 

In summary, although the 1/r 3 singularity makes 
rcnormalization process mandatory to implement the 
physical requirement of short distance insensitivity, the 
important aspect here is that this requirement remains 
equally valid even if there are no singularities at all. 
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